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ABSTRACT 

In this work, we introduce the use of H i Pfund (3 (Pf/3; 4.6538 /im) as a tracer of mass accretion from 
protoplanetary disks onto young stars. Pf/3 was serendipitously observed in NIRSPEC and CRIRES 
surveys of CO fundamental emission, amounting to a sample size of 120 young stars with detected Pf/3 
emission. Using a subsample of disks with previously measured accretion luminosities, we show that 
Pf/3 line luminosity is well correlated with accretion luminosity over a range of at least three orders of 
magnitude. We use this correlation to derive accretion luminosities for all 120 targets, 65 of which are 
previously unreported in the literature. The conversion from accretion luminosity to accretion rate 
is limited by the availability of stellar mass and radius measurements; nevertheless, we also report 
accretion rates for 67 targets, 16 previously unmeasured. Our large sample size and our ability to 
probe high extinction values allow for relatively unbiased comparisons between different types of disks. 
We find that the transitional disks in our sample have lower than average Pf/3 line luminosities, and 
thus accretion luminosities, at a marginally significant level. We also show that high Pf/3 equivalent 
width is a signature of transitional disks with high inner disk gas/dust ratios. In contrast, we find 
that disks with signatures of slow disk winds have Pf/3 luminosities comparable to those of other 
disks in our sample. Finally, we investigate accretion rates for stage I disks, including significantly 
embedded targets. We find that stage I and stage II disks have statistically indistinguishable Pf/3 line 
luminosities, implying similar accretion rates, and that the accretion rates of stage I disks are too low 
to be consistent with quiescent accretion. Our results are instead consistent with both observational 
and theoretical evidence that stage I objects experience episodic, rather than quiescent, accretion. 
Subject headings: stars: pre-main sequence — stars: planetary systems: protoplanetary disks matter 



1. INTRODUCTION 

In the study of protoplanetary disks and protostars, 
much effort has been focused on the study of mass accre- 
tion rates — the rates at which mass is transferred from 
circumstellar disks to stars — because it is so intricately 
linked to processes important for star and planet forma- 
tion. Mass accretion is a measure of the viscosity of the 
disk and determines the overall rate of mass and momen- 
tum transfer, and thus the pace of disk evolution. The 
rate of mass accretion will affect the disk lifetime (and 
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thus the time available for planet formation) and the 
rate of planetary migration, and may in turn be a tracer 
of th e presence of planets (e.g., [Alexander fc Armitagd 
l2007f l. Accretion greatly affects the inner disk envi- 
ronment, with the disk truncated and material lofted 
onto the star a t or near the stellar corotation radius 
(|Shu et al.lll994[) , which could be related to th e observed 
pile-up of giant planets at small orbital radii ()Lin et al.l 
119961 : iButler et ahl 120061 ). Finally, accretion is a tracer 
of star/disk magnetic interactions and determines the 
early angular momentum ev olution of the star (e.g., 
lAgapitou fc Papaloizou|[2000[ ). 

The most accurate measurement of mass accretion is 
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likely obtained from the spectrosco pic observation and 
modeling of UV exc e ss flux (e.g.. IValenti et al.l [T993I : 
IGuUbring et"all [19981: iHerczeg fc HiUenbrandl I2008ir as 
this provides a nearly direct measure of the total accre- 
tion luminosity. However, there has long been an interest 
in measuring accretion rates with data that are relatively 
simpler to obtain and analyze, as well as data that can 
be obtained at longer wavelengths, where extinction is 
lower. This has led to a number of studies of easily- 
observable H I emission lines believed to be produced in 
the accretion columns a nd accretion shock along w ith the 
UV continuum excess (jCalvet fc GullbrindlT998h . Such 
studies have shown that H i line luminosities correlate 
with UV-excess-derived accretion luminosities, and thus 
can be used as reasonably reliable tracers of mass accre- 
tion rates (e.g., Mu zcroUc ct al. 1998; Natta ct al. 2004 
Garcia Lopez et al.l 120061: INatta et al.l 120061 : IFang et al] 
200911 . albeit with some systematic unce rtainties and 
caveats (e.g.. IHerczeg fc HiUenbrandl I2008D . The ease of 
collecting spectra of H i emission lines results in more 
comprehensive samples of mass accretion rates, allows 
for the study of accretion rates in more embedded disks, 
and allows for a simultaneous measure of accretion rates 
and other disk properties, such as veiling by the disk con- 
tinuum, or molecular emission line strengths. In fact, one 
of the best studied H l lines — Ha — remains a popular 
means of estimating mass accretion, in spite of the fact 
that model s suggest it saturates e ven at moderate accre- 
tion rates (Muzer olle et al.l 11998"). and the fact that its 
strength and profile shape can also be affected by other 
para meters, including ou tflow rate and system inclina- 
tion ([Kurosawa et al 1120061 ). In addition, since measured 
accretion rates span several orders of magnitude, even 
measurements with error bars as large as ~ 0.5 — 1 dex 
can be used to broadly characterize a sample of objects 
and provide scientifically useful information. 

In this work, we introduce a new tracer of mass ac- 
cretion — specifically, H i Pf/3 — that offers several ad- 
vantages over other tracers, and results in one of the 
most comprehensive coherent datasets of a single mass 
accretion tracer to date. H i Pf/3 (n = 7 5; here- 
after Pf/3) is located at 4.6538 ^m in the M band, 
which places it between the R(0) and R(l) lines of 
the CO ro-vibrational fundamental band. Thanks to 
two large campaigns designed to study CO fundamental 
emission in protoplanctary disks with Kcck-NIRSPEC 
(e.g., Blake & B oogcrt 2004; Salvk ct al. 2 011) and VLT- 
CRIRES (e.g.. [Pon toDpida net al. l l2011at IHerczeg et all 
[20111 : iBrown et al.l Isubmittedl ) . Pf/3 has been serendipi- 
tously observed in more than 100 young stars with disks. 
Pf^ has similar equivalent widths to other infrared trac- 
ers, including Br7 (n=7— >4), with which it shares the 
same upper level energy, and appears to be ubiquitous 
for disks determined by other means to be typical ac- 
tive accretors. It also offers several distinct advantages 
over other accretion tracers due to its position in the 
M band. Firstly, it is an ideal tracer of accretion rates 
in transitional disks, as the low M-band photospheric 
flux results in a very high Pf/3 equivalent width when 
the dust continuum is low, as it is in transitional disks. 
Secondly, the high degree of veiling in the M band rel- 
ative to shorter wavelengths makes contamination from 
stellar photospheric absorption negligible in nearly all 
cases. Thirdly, tentative discoveries of young planets in 



disks ([Huelamo et al.ll2011l: iKraus fc Irelandll2012D have 
sparked interest in relating CO emission line variabil- 
ity to the tidal influence of proto-planets ([Regalv et al.l 
120111 ). An understanding of the connection between ac- 
cretion and Pf/3 may help observers differentiate between 
accretion- and planet-induced variability observed in CO 
emission lines. Finally, Pf/3 is observable even in heavily 
extincted disks, and so can be used to study the youngest 
disks embedded in their natal cloud. 

The ability to measure accretion rates in embedded 
disks is a particularly important strength of Pf/3 as an 
emission tracer. A star is thought to gain most of its mass 
while it is still embedded in a dense envelope, prompt- 
ing the naive expectation that accretion rates should 
be higher during these stages. During this phase, how- 
ever, the commonly-used accretion indicators in the op- 
tical and near-IR are not usually observable and this 
expectation is difficult to confirm. As a substitute for 
direct measurements of accretion rates, the mass his- 
tory of stars has been inferred fro m a global analysis 
of temperature - luminosity diagrams (iKenvon et al.lll990l : 
iDunham et aH [Mot iZhu et al.l l2010f ). In these anal- 
yses, the luminosity distribution is ~ 10% of the ex- 
pected luminosity distribution required to build stars 
via steady accretion, suggesting that the star may need 
to grow mostly in large, short bursts. A few studies of 
optical and near-IR accretion tracers in more extincted 
disks (.M uzcroUc ct al. 1998; White & Hillenbrand 200| 
[White et al.l 120071 ) have suggested a similar result — 
namely, that accretion rates in younger disks are lower 
than expected. However, as even near-IR accretion trac- 
ers like Ha, Br7 and Pa/3 are not visible in the most 
embedded disks, these studies may have been biased to- 
wards older systems, including possibly edge-on, evo- 
lutionarily older (yet observationally extincted) disks. 
Therefore, the study of Pf/3 provides an exciting new 
path for measuring accretion in the youngest, most em- 
bedded disks. 

In this paper, we report the detection and measure- 
ment of Pf/3 emission lines observed with Keck-NIRSPEC 
and VLT-CRIRES. We develop a method to use Pf^ 
emission to measure accretion luminosity by correlating 
line luminosities to known accretion luminosities for a 
sample of young stars. This correlation is then applied to 
a large sample of targets to provide accretion luminosi- 
ties for 120 young stars, including a sample of objects 
still embedded in their molecular envelopes. In Section 
U we briefly describe the NIRSPEC and CRIRES ob- 
servations and data reduction procedures, as well as the 
sample selection. In Section (U we discuss the flux ex- 
traction procedure and corresponding uncertainties. In 
Section Ism we calculate a relationship between Pf/3 line 
luminosity and accretion rate, and in Section [5.21 we ap- 
ply this to our full sample to provide accretion luminos- 
ity estimates for 120 stars. In SectionjBl we discuss some 
implications of our results. 

2. OBSERVATIONS AND REDUCTION 

2.1. Observing Procedures 

NIRSPEC observations were obtained with the Keck 
II telescope as part of a large survey of CO rovibra- 
tional emission from young stars with disks (see, e.g., 
ISalvk et al...201L) . spanning 2001-2011. Spectra were 
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obtained in the high- resolution mode with R~25000, in 
the M-wide filter. Echelle and cross-disperser positions 
were optimized for the observation of the P branch of 
CO v= 1^0. Although the favoured echelle settings 
evolved somewhat throughout the course of the survey, 
nearly all sources were observed in a setting that in- 
cluded the H I hues Pf;3 (4.6538 fj.m) and Hue (4.6725 
^m). NIRSPEC targets were observed in ABBA nod 
sets, with AB pairs differenced to remove thermal emis- 
sion from Earth's atmosphere. Integration times (expo- 
sure time multiplied by coadds) in each position were 
limited to 1 minute to minimize atmospheric variation 
between frames of a pair. To correct for telluric absorp- 
tion, A or B telluric standard stars with nearly blackbody 
spectra were observed close in time and airmass to the 
targets. 

CRIRES observations were obtained with the VLT as 
part of a large survey of CO rovibrat ional emission from 
young stars with circumstellar disks (jPontoppidan et aTl 
l2011at iBrown et al.l Isubmittedl ). i ncluding a significan t 
sample of embedded protostars (|Herczeg et al.l 120 111 ). 
Only ^0.02 /im are covered in each integration with the 
non-cross-dispersing CRIRES instrument, but the ma- 
jority of targets in this survey were still observed in a 
setting that included Pf/3. Observations were obtained 
with the 0".2 slit, resulting in a resolution of R '^90,000. 
Targets were observed in ABBA nod sets, with 60 sec- 
ond integration times for each image. Random dithers 
between and 1" were added to the standard 10" nod 
sequence to distribute the counts over different pixels, 
and ensure that the data would not always land on a 
bad pixel. When possible, brighter targets were observed 
with adaptive optics. 

A log of the NIRSPEC and CRIRES observations in- 
cluded in this study can be found in Table [TJ 

2.2. Data Reduction 

NIRSPE C data reduction routin es were developed by 
our team ()Blake fc BoogertI [200l in IDL. NIRSPEC 
spectra were first linearized using the shape of telluric 
emission lines and the observed spectra. Fluxes were ex- 
tracted with a 2.8cr aperture around the best-fit PSF cen- 
ter, and nearby columns were subtracted to account for 
any residual sky emission in the A-B difference image. 
Extracted spectra were wavelength calibrated utilizing 
telluric lines. Telluric standards were processed with the 
same procedures as the target stars. Target spectra were 
then divided by the telluric standard star spectra (ad- 
justed for small differences in airmass) to obtain spectra 
corrected for telluric absorption. 

Although telluric standards were nearly featureless, 
most A stars have photospheric H i absorption, including 
that of Pf/3. Thus, correcting for this feature is crucial for 
obtaining accurate Pf/3 fluxes in our targets. To account 
for this, telluric standard spectra were first flattened uti- 
lizing Kurucz models, broadened to match observed H i 
absorption shapes. An example of this process is shown 
in Figure [H Note that the corrected telluric spectrum 
should be linear (except for the narrow telluric absorp- 
tion lines), and that this can easily be inspected by eye to 
confirm proper correction of the photospheric H I. In rare 
cases, the H i absorption profiles were not well accounted 
for by Kurucz models; in these cases, we simply allowed 
for an additional Gaussian-shaped absorption component 
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Figure 1. Solid and dotted lines show sample standard-star spec- 
trum before and after stellar atmospheric model division. The 
shaded region marks an approximate location for Pf/3 emission, as- 
suming a width of 190 km s~^. The actual region will depend on 
the line width and relative radial velocity between the source and 
the standard stars. 

to be removed from the spectrum. For all spectra, we 
performed an additional check on potential contamina- 
tion from the standard star: we divided a flat standard 
star observation by each Kurucz-modcl-correctcd stan- 
dard star observation to confirm that it did not produce 
a spurious H i emission line. 

Since the primary observational target of these sur- 
veys, CO, can overlie telluric CO absorption lines, tar- 
gets were typically observed on two or more dates, and 
spectra combined to create complete CO line profiles. 
The H I line emitting regions are relatively unaffected 
by telluric absorption and so the repeated observations 
simply increase the S/N in these lines. Alternatively, 
multiple epochs can also be left separate to investigate 
accretion variability. This aspect of the dataset is left 
as future work, and will likely benefit from the inclusion 
of additional datasets from the NIRSPEC and CRIRES 
archives. 

Approximate fiux calibration was achieved by compar- 
ing telluric standard star spectra with literature photom- 
etry to derive a conversion from counts to flux. In all 
cases, the brightest observation of standard and target 
were utilized, as these would represent the best-centered 
observations. However, there is no guarantee that the 
sources were well centered, and the flux correction factor 
due to slit losses can be either greater than or less than 
one, depending on whether it is the source or the telluric 
standard that is observed off center. Nevertheless, we 
find decent agreement between absolute flux and litera- 
ture photometry; the RMS difference between measured 
and literature fluxes (see Table [5]) is 50% for our sample 
as a whole. 

Due to the inherent difficulties with absolute flux cali- 
bration for spectroscopic observations, we have chosen 
to correct observed fluxes with literature photometry 
whenever possible. The photometry used in this work 
can be found in Table [H We used M-band photom- 
etry when available, or secondarily an interpolation of 
other measurements, often Spitzer IRAC 4.5 and 5.8 /xm 
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fluxes. Although young stars with disks are known to 
have variable infrared emission, the vari ability is typi- 
cally less than 10% for class O-II disks ()Luhman et alj 
[2nTnh . Va riability may be high er for transitional disks 
(^30-50%: lEspailjat et a r 2011^ and outbursting sources 
like EX Lup fe.g. JAspin e t al. 2010), and so we note that 
such sources can have a corresponding systematic error 
in their Pf/3 line flux measurement. If photometric mea- 
surements were not available in the literature, we sim- 
ply used our own absolute photometry derived from the 
spectroscopic observations. 

CRIRES reductions util ize routines written by K. Pon- 
toppidan and described in lPontoppidan et al.l ()2008D . 2D 
images for each source were added together after aligning 
and resampling onto a 2x finer grid. Spectra were ex- 
tracted using an aperture twice the size of the FWHM of 
the spectral profile. As with the NIRSPEC observations, 
nearby telluric standards were observed close in time and 
airmass. Wavelength calibration was performed on the 
telluric standards using telluric emission lines and was 
applied to the target spectra after small shifts to align 
the spectra of target and standard. 

In the nominal CRIRES reduction routines, target 
spectra were then simply divided by the telluric spec- 
tra to remove telluric absorption features. However, this 
does not account for H l absorption in the standard stars, 
and would have produced anomalously strong H i emis- 
sion in our targets. Therefore, we implemented a routine 
adapted from the NIRSPEC reduction routines, which 
uses Kurucz photospheric models to fit and correct for H 
I absorption features in the standard stars. 

Spectra of standard stars were also compared with lit- 
erature photometry to obtain a conversion from counts 
to pixel, which were applied to the target spectra to ob- 
tain an approximate absolute flux calibration. However, 
just as with the NIRSPEC data, CRIRES flux calibra- 
tion suffers from an uncertain degree of slit losses for any 
given observation. Therefore, as with the NIRSPEC ob- 
servations, we adjusted the spectra to match literature 
photometry whenever possible. We flnd a few large dis- 
crepancies (factors of ^ 2 — 3 between literature and mea- 
sured photometry, including for Elias 29, IRS 44, LLN 
8, LLN 17 and WL 12. Since all show higher litera- 
ture fluxes than CRIRES-measured fluxes, it is possible 
that this is because prior observations with the Spitzer 
Space Telescope or ISO included more than one source in 
their relatively large spatial point-spread-function. For 
the non-outliers, we flnd the RMS difference between ob- 
served and literature fluxes to be ^52%. 

2.3. Sample Selection 

The NIRSPEC and CRIRES target lists used in this 
work were compiled for prior studies of M-band CO emis- 
sion and do not represent an unbiased sample of young 
stars. However, the large sizes of these surveys means 
they come close to representing a complete flux-limited 
sample of disks with detectable CO fundamental emis- 
sion in nearby clouds. The NIRSPEC M-band survey as 
a whole is dominated by revealed, optically thick disks, 
including significant numbers of both low-mass (T Tauri) 
and mid-mass (Herbig Ae/Be) disks. It also includes a 
smaller but significant number of embedded disks. The 
NIRSPEC sample is strongly biased against sources with 
tenuous disks or no disks at all, as these disks tended not 



to produce observable CO rovibrational emission lines. 
The NIRSPEC sample is also biased towards clouds at 
the higher declinations observable from Mauna Kea, sam- 
pling well in Taurus, Perseus, Serpens and Ophiuchus, 
only poorly in Lupus, and not at all at lower declina- 
tions. The NIRSPEC sample was also limited to targets 
with M-band continuum fluxes brighter than ^0.01 Jy 
(M-9). 

The CRIRES M-band CO survey contains many of the 
bright Class II disks visible from the southern sky as well 
as a number of Class I embedded protostars and transi- 
tional disks. This includes targets in Lupus, Vela and 
Chamaeleon that are difficult or impossible to observe 
with NIRSPEC. 

A more detailed discussion of the evolutionary status 
of our sample as a whole is included in Section 16.41 

3. GENERAL CHARACTERISTICS OF PF^ DETECTIONS 
AND NON-DETECTIONS 

A sample of targets showing different strengths and 
types of Pf/3 emission profiles is shown for NIRSPEC in 
Figure H and for CRIRES in Figure E] Although the 
biases in our sample selection make any statistics diffi- 
cult to interpret, Pf/3 appears to be a robust tracer of 
accreting systems. Greater than 80% of optically thick, 
classical disks around T Tauri and Herbig Ae/Be stars 
show detectable Pf/3 emission. It is interesting to note 
that we find a slightly lower detection fraction for embed- 
ded protostars 60%). A similar fraction of embedded 
sources without Pf /3 emission is found in a lower resolu- 
tion VLT-ISAAC study of ^30 sources by Pontoppidan 
et al. (2003, their Figures 1-7). Although the differ- 
ence in detection rates may not be significant, especially 
because of possible biases in the selection of embedded 
objects, the presence of a large number of embedded pro- 
tostars with little or no detectable Pf/3 is interesting in 
and of itself. Because embedded protostars have rising 
continua towards the infrared, one might imagine that 
this could decrease line/continuum ratios and introduce 
an obsevational bias. However, the embedded protostars 
with Pf/3 detections do not appear to be systematically 
fainter than those without, nor do the embedded proto- 
stars appear to have systematically higher M-band fluxes 
than the protstars with classical disks. 

The majority of targets with Pf/3 non-detections rep- 
resent a few distinct types of objects. Not surpris- 
ingly, debris disks and disks classified as weak-line T 
Tauri stars (wTTS's) or protostars with class HI spec- 
tral energy distributions (SEDs) generally do not appear 
to show Pf/3 in emission. wTTS's with weak or non- 
existent Pf/3 emission include HD98800 (also a circumbi- 
nary disk; Furlan et al. 2006), Hen 3-600A, LkHa 332 
Gl, TWA 7, TWA 8a and WaOph 4. These additionally 
show no CO fundamental emission and very low levels 
of veiling, likely reflecting the expected relationship be- 
tween warm inner disk gas, inner disk dust, and accre- 
tion. Other targets without detectable Pf/3 emission are 
in some way straddling the evolutionary boundary be- 
tween optically thick and tenuous disks. These include 
HD 36917 (classified as trans itioning between classes II 
and HI inlManoi et all [2001 . CoK u Tau/3 (class II i n 
[Andrews k Winiamsll2005r wTTS in lFurlan et alllMll) . 
FN Tan (classified as bordering between classical T Tauri 
star — cTTS — and wTTS in .Furlan et al..,2011, ). HQ 
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Figure 2. Sample set of Pf/3 emission line profiles from full NIR- 
SPEC sample. (Note that any overlying CO emission or absorption 
lines have already been removed.) 

Tau (wTTS by Ha equivalent wid th definitions, but 
cTTS according to its Ha 10% width iFurlan et aLll201lD 
and WSB 60 (classified as Class H in lEvans et al.l 120091 
bu t shows transit i onal disk millimeter emission shape 
in [Andrews et al.l 1201 Other targets with no de- 
tectable Pf/3 emission include the transitional disks HD 
149914 (which also shows no detec table Br7 eniission ; 
iBrittain et "alll20 07i as well as SR 21 (jBrown et al.ll2007D . 
SR 21 has detectable CO fundamental emission, although 
in contrast to many other transitional disks, the emis- 
sion arises from moderately large disk radii (^ 5-8 AU; 
iSalvk et al.|[20Tlt iPontoppidan et all 1201 Ibl) . Two cir- 
cumbinary disks, CoKu Tau/ 4 (7.8 AU ; Irela nd fc KrausI 
[2008h and ROXs42c (23 AU: IKraus et al.ii20nD . also do 
not show Pf/3 emission, although it should be noted that 
this docs not appear to be universally true, as close bi- 
nary DP Tau (projected separation 15.5 AU IKraus et al.l 
l201lb has detectable Pf^ emission. A final interesting 
set of objects that universally shows no detectable Pf/3 
emission is the set of FU Orionis stars. 

It is interesting to ask whether any classical, opti- 
cally thick disks show no detectable Pf/3 emission, ei- 
ther to test whether Pf/3 could be unreliable as an ac- 
cretion tracer, or to search for unique targets with op- 
tically thick disks but no accretion. The list of disks 
with class II SEDs that do not show Pf/3 emission above 
detectable limits include c2dJ033035.92-h303024.4, HK 



Figure 3. Sample set of profiles from full CRIRES sample. (Note 
that any overlying CO emission lines have already been removed.) 

Tau, IRAS 03380-^3135, IRAS 04385-1-2550, LkHa 270, 
LkCa 8, LkHa 325 and WaOph 5. The first five show no 
apparent CO emission, while the final three targets show 
weak CO fundamental emission, and their spectra sug- 
gest the possible presence of Pf/3 at low S/N; therefore, 
these suggest the presence of inner disk gas. Few disks 
are seen with strong CO fundamental emission and no 
detectable Pf/3 emission, one exception being the tran- 
sitional disk SR 21. Very low accretion rates are also 
apparently not detectable in our sample. For example, 
the spectroscopic binary Hen 3-600A, accretiii g at a rate 
of ~ 5 X 10-"MQyr-i ijMuzeroUe et al.|[2000l ). does not 
have detectable Pf/3 emission. 

We investigated but did not determine any relationship 
between line shapes or widths and any other stellar or 
disk parameters, except that sources with double-peaked 
fine profiles (including CI Tau, HD 141569 A, RY Tau 
and SU Aur) tend to have moderately high disk inclina- 
tions (typically 50-60°). A full discussion of the relation- 
ship between CO vibrational emission and Pf/3 emission 
is left as future work. Here we focus on targets with de- 
tectable Pf/3 emission, in order to measure mass accretion 
rates. Therefore, for this study we have selected a sub- 
set of the full NIRSPEC and CRIRES samples via visual 
confirmation of the presence of Pf/3 emission. The com- 
plete target list used in this study, with relevant stellar 
parameters, is shown in Table |3l 
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Pf/3 lines are unresolved in all natural-seeing NIRSPEC 
observations, constraining the emission to radii less than 
~50 AU for typical seeing and stellar distances. Pf/3 
observations are also unresolved in the adaptive optics 
corrected CRIRES spectra. With a typical point spread 
function core of 0.18", this constrains the Pf/3 to radii 
less than ~ 13 AU at 140 pc. This is consistent with an 
accretion origin for Pf/3, as opposed to an outflow origin. 

4. DERIVATION OF PF/3 LINE LUMINOSITIES 

4.1. CO Emission Corrections 

The Pf/3 line center is close to a number of CO vibra- 
tional lines, including v — 1 ^ Q R(0) (4.6575 /im) and 
R(l) (4.6493 ^xm), v = 2-^l R(8) (4.6523 ^m) and "CO 
R(14) (4.6572 ^m) and R(15) (4.6504 /xm). A majority 
of Pf/3 profiles are affected by at least one of these CO 
features, which contain non-negligible amounts of flux. 
A number of techniques were utilized to correct for the 
CO emission, as demonstrated in Figure SI If possible, 
CO lines were removed via subtraction of a Gaussian fit 
to the CO emission lines. In cases where the fits were 
not reliable, the line was replaced by a linear interpola- 
tion of the surrounding flux. In some cases, the nearby 
w = 2 — >■ 1 or ^•^CO lines contributed so significantly to 
the flux that it was difficult to separate the CO and Pf/3 
contributions. In such cases, an adjacent CO line would 
be shifted to the contaminating line location and sub- 
tracted. This necessarily introduces some error (of order 
10%) due to the fact that the flux in adjacent CO lines 
may not be identical. However, more precise modeling 
to derive these CO line fluxes was not pursued, due to 
severe blending of the v — 2 1 and ^'^CO rotational 
ladders with the much stronger v — 1 Q emission lines. 

4.2. Flux Extraction 

Observed Pf/3 line profiles are typically not described 
by a Gaussian or other simple functional form, and so 
we chose to extract fiuxes by summing all fiux within 
a defined window, after subtraction of the continuum. 
The extraction window and continuum level were both 
determined by eye, and have associated systematic un- 
certainties. These can be particularly high if a number of 
adjacent CO lines makes the continuum hard to define, 
and/or if the Pf/3 profile has broad wings. The sam- 
ple overlap between NIRSPEC and CRIRES allows us 
to estimate this uncertainty to some degree, by compar- 
ing the fluxes derived from each dataset. We flnd good 
agreement between extractions using the two datasets, 
with typical differences less than a factor of two. Note 
that since the Pf/3 line luminosity or the underlying con- 
tinuum can be variable, some of this difference could be 
attributable to real changes that occurred between the 
two sets of observations. In our most egregious case, HL 
Tau, we flnd a factor of 3 discrepancy between the two 
datasets. This extreme case is illustrated in Figure [5l 
showing how uncertainty in both the continuum position 
and the position of the Pf/3 line wings results in large 
flux differences between the two extractions. 

A few of the photometrically-corrected NIRSPEC line 
proflles are shown in Figure IHl the remaining NIRSPEC 
and CRIRES proflles are available in the on-line version 
of the article. Our derived Pf/3 luminosities, computed 
using the distances in Table [3j are listed in Table H] Line 



luminosities have also been correct ed for extinction, as- 
summg an extinction of 0.034 x Ay (iCardelh et al.lll989D 
at 4.6538 /xm, and the visual extinctions in Table |3l if 
the extinction is not known, we do not include any cor- 
rection. Note that even the maximum extinction in our 
sample {Ay — 34) only results in a change in line flux 
by a factor of 3, and thus computation of Pf/3 lumi- 
nosity is relatively insensitive to extinction. When both 
NIRSPEC and CRIRES data were available, we show the 
average of those measurements. 

4.3. Importance of Photo spheric H I Absorption 

While accretion flows are believed to be the source of 
the Pf/3 emission, the observed spectra also include a con- 
tribution from the underlying stellar photosphere of the 
star/disk system. As we will discuss in this section, the 
correction for the underlying stellar photosphere should 
in nearly all cases be insigniflcant compared to other un- 
certainties, and we therefore do not apply any correction 
to most sources. (Note that this correction is distinct 
from the correction for Pf/3 absorption in the photosphere 
of the telluric standard star, which is always corrected for 
in the data reduction process). 

Stellar photospheres have H i absorption features that 
will tend to reduce the amount of observed Pf/3 flux, i.e. 
Pf/3 fluxes must be increased somewhat to correct for 
stellar absorption. The degree of correction depends sen- 
sitively on the stellar spectral type and veiling. Figure 
[71 panels a-e, show synthetic spectra and observed spec- 
tra for five spectral types (AOV ~ MOV), demonstrating 
that the absorption is most prominent at early spectral 
types. However, M-band veiling values for optically thick 
disks act in the opposite sense, with much higher veil- 
ing around early-type stars. Thus optically thick disks 
have negligible corrections to Pf/3 from underlying stellar 
absorption. Synthetic veiled spectra with veiling values 
typical of optically thick disks (r = 46, 18, 10, 7 and 3 for 
spectral types A-M, derived from disk SED s produced 
with RADMC; iDullemond fc DominikI |2004[ ) are shown 
as dotted lines, and are indistinguishable from straight 
lines. 

Stars with G~M spectral types have M-band spectra 
that are dominated by a molecular pseudo-continuum, 
rather than simply H i absorption (see panels c, d and 
e), which can result in some change to Pf/3 flux due to 
coincident molecular absorption features near Pf/3. The 
equivalent width of the underlying photospheric features 
drop by a factor of 1 -f r where r is the veiling (deflned 
as the dust continuum flux density divided by the stellar 
photospheric flux density). Considering the spectrum 
of HD 79210 shown in Figure [71 a reasonable M-band 
veiling value of ^3 for an optically thick disk reduces the 
strongest photospheric features to the 7% level; this 
and their narrow width results in a negligible correction 
for optically thick disks. 

Transitional disks — disks with inner regions depleted 
of small dust grains — have low continuum fluxes, and so 
could potentially have large relative contributions from 
the underlying photosphere. However, we flnd that such 
targets (see Figure [3 panels f), h) and j)) often have high 
line/continuum ratios, likely because in spite of their low- 
ered near-IR continuum fluxes, tra nsitional disks hav e 
only slightly lowered accretion rates (iNaiita et al.ll2007l) . 
(See extended discussion in Section [6]). As examples. 
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Figure 4. Gray and black lines show sample spectra before and after removal of CO emission lines, respectively. Dotted lines mark CO 
1—^0 lines; dashed lines mark CO 2—^1 lines, (a) Lines are removed with Gaussian fits, (b) R(0) is removed with a linear fit to the 
underlying spectrum, (c) R(8) is removed by using a Gaussian fit to R(7). 

molecular absorption features in its stellar photosphere 
have a non-trivial effect on the determination of the disk 
Pf/3 line flux. We describe this example in detail in an 
Appendix. 

The largest possible correction would be required for a 
naked A or F star with a low Pf/3 line/continuum ratio 
(see, e.g., panel g), perhaps doubling the ratio of true to 
observed Pf/3 emission flux. Such a target would be diffi- 
cult to identify because its relatively featureless M-band 
spectrum looks the same whether it is naked or veiled. 
To our knowledge, there are no known diskless targets in 
our sample that also show detectable Pf/3 emission. Of 
course, stellar Pf/3 absorption could erase signatures of 
weak Pf/3 emission from weakly- accreting disks; however, 
we have no good way to pinpoint such targets. A possible 
example of such a target could be SR 21 (see Figure [5]) a 
G2.5 star with a transitional disk, which has CO gas at 
~7AU but no measurable accretion (jPontoppidan et aD 
[20081) . 




4.645 4.650 4.655 4.660 

Wavelength [yam] 

Figure 5. Pf/3 line profiles from HL Tau observed by NIRSPEC 
and CRIRES before (black) and after (gray) CO emission removal. 
Dotted lines mark the regions chosen for flux extraction, demon- 
strating how uncertainty in the continuum level and line width can 
in extreme cases result in a factor of three difference in extracted 
line flux. 

corrections to the Pf/3 flux due to photospheric H i from 
transitional disks HD 141569 A, LkHa 330 and GM Aur, 
(panels f, h and j) would be only ^ 20, 45 and 3%, re- 
spectively, even in the worst case scenario (no veiling). 
Actual corrections for these sources should be a factor of 
a few lower, due to non-zero M-band veiling fSal vk et al.l 
I2009D . The molecular photospheric absorption features 
might be more problematic, but luckily the spectra can 
be examined at other wavelengths to determine whether 
there is likely significant contamination from the under- 
lying photosphere. With visual examination, we found 
one target — DoAr 21 — which shows both low Pf/3 
equivalent width and low veiling, and in this case the 



5. PF/3 LUMINOSITY AND ACCRETION RATES 

5.1. Correlation between Pfj3 Luminosity and Accretion 
Luminosity 

In Figure [U we show the correlation between our mea- 
sured Pf/3 line luminosities and previously-measured ac- 
cretion luminosities, listed in Table [5j We chose to 
maximize the number of sources in this correlation and, 
therefore, the literature accretion luminosities are de- 
rived from a variety of observational methods. The ma- 
jority are derived from models and observations of the 
UV excess continuum, while others are derived some- 
what less directly from UV data, utilizing a relation- 
ship b etween mass accretion and the Balmer disconti- 
nuity (jMuzerolle et al.ll2004[ ) or approximating the UV 
excess using photometry. About a third of the accretion 
rates are derived from empirical relationships between H 
I emission lines (Br7 or Pa/3) and accretion luminosity. 
Although the Br7 and Pa/3-derived accretion rates are 
themselves an indirect measure of accretion rate, we find 
no significant biases between these and the UV-derived 
accretion luminosities in Figure [51 

Because the accretion rates of high-inclination sources 
can be unreliable, we remove these from our correlation. 



8 



Salyk et al. 




Figure 6. Pf/3 line profiles observed with NIRSPEC, after correction for CO emission and literature photometry, when available. Dotted 
lines mark the baseline and limits used for flux extraction. The full set of NIRSPEC and CRIRES line profiles is available in an online 
version of this figure. 
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Figure 7. Panels a)-e) show synthetic stellar H I absorption spectra without (solid curves) and with (dotted curves) the addition of typical 
veiling from an optically thick disk. Q Panels a), c), d) and e) additionally show observed spectra of unveiled stars, with the K4 and MO 
spectra demonstrating that the spectra of cool stars are dominated by a molecular pseudo-continuum. Panels f)-j) show example spectra 
(black) and worst-case Pf/3 flux corrections (i.e., assuming no veiling; gray) for the following targets: f) transitional disk HD 141569 A, g) 
Herbig Ae/Be star HD 144432 S, h) transitional disk LkHa 330, i) cTTS DoAr 24 ES and j) transitional disk GM Aur. 

''AOV and FOV spectra arc produced from ATLAS models using the ATLAS, WIDTH and SYNTHE Linux port fKuruc^ IT99^ : 
[Sbordone et al. 2004, 2005). GOV, KOV and MOV models are produced with the MOOG stellar synthesis code (Snedon 1973) using 
MARCS atmospheric models (iGustafsson et al]|2008l ). but show only the contribution from H I lines. The actual spectra are dominated 
by molecular absorption, and stellar models in this wavelength range have not been w ell benchmarked. All spectra assume a rotational 
broadening of 10 km s~^. Representative veiling values were derived using RADMC l IDuUemond Ik Dominiki i2004 ) . assuming optically 
thick disks that extend to Tin = 1500 K. 
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Figure 8. Correlation between literature accretion luminosity (see 
Tal3le[3ll and Pf/3 luminosity (this work). If Pf/3 is measured by 
both NIRSPEC and CRIRES, we show two symbols (squares and 
triangles, respectively) on the plot. Sources with high inclination 
or otherwise unreliable accretion rates (in gray) are excluded from 
the correlation analysis. The solid line shows the best linear fit to 
the data excluding points from Mendiguti'a et al. (2011) (Equation 
[TJ . The dashed line shows the best linear fit including points where 
La.cc is derived from Mondigutfa et al. (2011). 



We also omit HL Tau from the correlation; although 
its disk appears t o be viewed at moderate inclination 
(|Kwon et al.ll201ll ). its high level of veili ng makes litera- 
ture accretion rates very uncertain fe.g.. lWhite fc Ghe3 
[2OOII) . 

We find that the correlation between Pf/3 line lu- 
minosity and accretion luminosity is heavily influenced 
by targets with h i gh ac cr etion luminosities derive d by 
iMendigutia et all (|2011| ). iMendiguti'a et alj ()201lD de- 
rived accretion rates in 38 HAeBe stars by estimating 
the Balmer excess using U- and B-band photometry com- 
pared with accretion shock models. Since accretion rates 
correlate with stellar mass, these accretion rates popu- 
late the upper right corner of Figure HI If these points 
are included in the fit, we find that our correlation then 
poorly reproduces previo usly-measured accretion rates 
for low-mass stars. Indeed. iMendigutfa et al.l()2011[ ) note 
a change in the slope of the correlation between accre- 
tion rate and stellar mass between low-mass and high- 
mass stars, and find larger accreti on rates than those 
determined from Br7 emission lines (jGarcia Lopez et al.l 
I2OO6D , implying that correlations between accretion trac- 
ers and accretion rates may not extend linearly to higher 
mass stars. However, UV excesses are significantly more 
difficult to measure around high-mass stars, and addi- 
tional study is warranted to understand whether the 
emission line-accretion luminosity relationship needs to 
be modified at higher stellar masses. Therefore, in this 
study, we omit the acc retion luminosities derived by 
IMendi gutia ct al. (2011) from our analysis. If the re- 
sults of Mcndigutfa ct al. (2011) are correct, our study 
will tend to underestimate accretion luminosities for 
> SMq. 

Excluding the accretion luminosities in 
IMendigutia et al.l ()201lD . we find the following relation- 
ship between Pf/3 luminosity and previously-measured 



Figure 9. Mass accretion rates derived from this work compared 
with values from the literature. High-inclination disks are shown 
plotted in gray. The dashed line shows a 1 : 1 correlation and dotted 
lines show a difference of 0.5 dex. 

accretion luminosities: 

logLacc[i0] = (O.91±O.16)xlogLpf/3[L0]-f(3.29±O.67). 

(1) 

5.2. New Accretion Rates 

We use Equation [1] to compute accretion luminosities 
for our entire sample. These are listed in Table |H An 
estimate of M is derived using the standard relationship 
between Lace and M: 

Lace = 0.8 GM,M/i?, (2) 

from lGullbring et al.l (flOOl . Af* and are taken from 
Table [S] If not shown in Table [3l i?* is derived from L^ 
and Ti, assuming L^, = AirR^aT^. We do not compute M 
if these parameters are not available. However, M could 
easily be computed at a later time when these parameters 
are measured. 

In Figure [SI we show accretion rates derived here com- 
pared to existing values from the literature (Table [SJ. 
The accretion rates derived here show no systematic bias, 
and most accretion rates are within 0.5 dex of previ- 
ous measurements. The standard deviation of the cur- 
rent and previous measurements is 0.77 dex. These re- 
sults are similar to those for other emission line accre- 
tion tracers, albeit over a somew hat more limited range 
of acc retion rates. For example, iHerczeg &: Hillenbrandl 
(|2008D quote standard deviations of 1.0, 0.71 and 1.1 dex 
between UV continuum excess and Ca 11 A8542, Ca 11 
A8662 and He I A8576, probing accretion rates as low 
as 10^ ^^M(7)yT~^ . In addition, IHerczeg fc Hillenbrandl 
(|2008[) estimate a factor of 4 (0.6 dex) random error in 
UV excess measures alone, due to errors in extinction 
and dista nce. Also, accretion r ates are known to be vari- 
able, with INguven et all (|200 9^ finding typical accretion 
rate variations of 0.35 dex, but variations commonly as 
high as 0.5 dex. Therefore, Pf/3 appears to be an ac- 
cretion tracer on par with existing tracers, with scatter 
not much larger than what would be expected from ac- 
cretion variability and errors related to the derivation of 
UV-based accretion luminosities. 
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5.3. About Upper Limits 

Computing upper limits on the Pf/3 Une flux, and thus 
accretion rate, for any given source is not straightfor- 
ward, and because these are not crucial for the work 
presented here, we do not make an attempt to compute 
them. Complications include the fact that the Pf/3 line 
shapes are complex, and their widths varied, so the cor- 
rect choice of assumed lineshape and width is not obvi- 
ous. In addition, sources with no detectable Pf/? emission 
are often less-heavily veiled sources, whose underlying 
photospheres must be modeled to derive an accurate Pf/? 
flux. We encourage those interested in an upper limit for 
a particular source in the NIRSPEC or CRIRES archives 
to contact the authors of this work to discuss whether 
reasonable assumptions can be made to derive limits for 
this source. 

Nevertheless, we can make some general observations 
about how sensitive Pf/3 measurements are to accretion 
rates. Assuming the Pf/3 line is a Gaussian with (Jwidth — 
0.001/im (derived from a flt to AB Aur's Pf/? emission 
line) , then with a continuum flux density of i^^ = 1 Jy, a 
signal-to-noise ratio (SNR) of 50, a Scr line peak, and for 
a 1 solar mass, 1 solar radius star, the minimum measur- 
able accretion rate is ~ 4 x lO~^M0yr~^. The minimum 
measurable accretion rate then scales roughly propor- 
tional to . For the integration times used for the 
acquisition of the data used in this work, Pf/3-derived ac- 
cretion rates do not appear to be as sensitive as p revious 
studies using Br7 or Pa/3 (e.g., EattaeTaDIlQOa). 

The lowest detected accretion rates in our sample are 
for the transitional disks TW Hya and DM Tau, with 
rates of 2.6 x 10"^ and 4.3 x lO~^M0yr~^, respectively. 
Pf/3 is especially well-suited for measuring small accre- 
tion rates from transitional disks, as the dust-cleared in- 
ner holes decrease Fi, and increase the line/continuum 
ratio (see Section 15^ for more detail). However, as the 
SNR per unit exposure time also decreases as de- 
creases, sufficient time must be expended to reach a rea- 
sonable SNR. 

6. ANALYSIS AND DISCUSSION 

6.1. Utility of PfP and Comparison to Other Accretion 

Tracers 

In this work, we have introduced the use of H i Pf/3 to 
measure mass accretion from protoplanetary disks. Al- 
though numerous accretion tracers have been developed 
prior to this work, the use of Pf/3 offers several advan- 
tages, which we outline here. 

One major advantage of using Pf/3 as an accretion 
tracer is that it comes "for free" and contemporane- 
ous with observations of CO fundamental emission. The 
NIRSPEC and CRIRES surveys from which we extract 
the data in this work provide a self-consistent sample of 
accretion luminosity estimates for 120 targets, and the 
NIRSPEC and CRIRES archives likely include many ad- 
ditional protoplanetary disk targets for which accretion 
luminosities could be estimated. To our knowledge, this 
is the first work to make extensive use of M-band spec- 
tra of protostars from both NIRSPEC and CRIRES — 
data which are all currently available in their respec- 
tive archives. The varied and large set of targets re- 
sulting from this combined dataset is a great demon- 
stration of the utility of these large archives. As future 



work, we plan to extend our analysis to investigate long- 
term (~year timescales) accretion variability for different 
classes of disks. 

The simultaneous observation of CO and Pf/3 also al- 
lows for a comparison between the accretion fiow and 
molecular disk gas at distances of ^ 0.1-1 AU from th e 
star (e.g., iPontoppidan et al.|[20TIbt iSalvk et al.ll201lD . 
which is traced by CO fundamental emission. The re- 
lationship between t hese two disk componen ts is being 
investigated, e.g., bv lBrown et al.l (jsubmittedl ) . The con- 
temporaneous nature of these measurements also allows 
for simultaneous measurements of disk gas and accre- 
tion rate as a function of time, although it should be 
cautioned that there can be changes in emission line 
fiuxes without correspo nding changes in accretion rate 
(e.g., IGahm et al.l 120081 ). The tidal infiuence of proto- 
planets is predicted to cause CO emission line variability 
(IRegalv et al.ll201ll) . but understanding the influence of 
proto-planets requires ruling out other possible sources 
of line variability, such as changes in accretion rate or 
accretion geometry. 

Pf/3 (n=7— >-5) is in many ways similar to Br7 (n=7— >4; 
2.160 /im), originating from the same upper level energy 
and appearing in the infrared. In Figure 1101 we compare 
Pf/3 EW with literature values of EW(Br7) as well as 
EW(Pa/3) (see Table 15). Typical values of Pf/3 equivalent 
widths are between 1/3 and 3 x EW(Br7) or EW(Pa/3), 
with EW(Pf/3) typically being simflar to EW(Br7) but 
somewhat smaller than EW(Pa/3). The lack of a strong 
correlation between tracers likely reflects the varied tem- 
peratures of the accretion column, the underlying pho- 
tosphe re , and the co ntinuum veili ng. [M uzerolle et al] 
' 1998|) . iCalvet et ahl (|2004) and iDoneh ew & Brittair 



201 If ) all find correlations between accretion luminos- 



ity and Br7 luminosity very similar to our Equation 
J Equation [T corn bined with the relationships in 
Muzerolle et al.lllSQSD and IDonehew fc BrittainI (|201lD 
predict somewhat h igher EW(Pfj3 ) than EW(Br7) , while 
the relationship in ICalvet et al.l ()2004D predicts some- 
what lower EW(Pf/3) as compared to EW(Br7) (with 
differences less than ^ 0.2 dex). Lineshapes of the two 
accretion tracers appear to be quite similar — typi- 
cally single-peaked with li ne widths near ^ 200 km s^^ 
(jCarcia Lopez et al.|[2006l ). HD 141569 A, on the other 
hand, is an example of a disk with both double-peaked 
Br7 and Pf/3 emission lines. Interestingly, however, sev- 
eral targets for which Pf/3 is clearly in emission show no 
Br7, or Br7 in absorption — for example, HD 142527, 
HD 142666, T CrA and TY CrA. Thus, Pf^ may be a 
more robust accretion tracer in targets with low near-IR 
veiling. 

Although equivalent widths are similar, Pf/3 offers 
some distinct advantages over Br7 and Pa/3. One ad- 
vantage is the relatively higher continuum veiling in the 
M band as compared to shorter wavelengths. Since the 
photospheric emission is a smaller fraction of the total 
flux at 5 /im than it is at shorter wavelengths, Pf/3 suf- 
fers less contamination and uncertainty from the under- 
lying photospheric H i absorption (as discussed in de- 
tail in Section 14. 3p . Being at longer wavelengths also 
makes Pf/3 less sensitive to extin ction. Using the redden- 
ing law of lCardein et al.l (| 19891) . A^i.^/Av = 0.12 while 
Apfp/Ay = 0.03, a factor of 4 difference in magnitudes 
of correction. At Ay ~ 30, for example, this results in 
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Figure 10. EW of Pf/3 compared to those of Br7 (bottom) and 
Pa/3 (top). The outliers with highest Pf/3 EW in the bottom plot 
are transitional disks HD 141569 A, GM Aur and TW Hya. 

a flux correction factor of ^ 28 for Br7 but only for 
Pf/3. Finally, Pf/3 equivalent widths are quite high for ac- 
creting disks with low continuum veiling; thus, Pf/3 is an 
excellent tracer for low accretion rates in disks with inner 
regions de plete d in small dust grains (discussed further 
in Section [6.2^ . One disadvantage of Pf/3 as an accre- 
tion tracer is the systematic errors that are introduced 
by complex CO emission/absorption spectra, which can 
result in up to a factor of a few uncertainty in flux in ex- 
treme cases (as in Figure [5]). Another is that the veiling 
is difficult to determine empirically in earlier-type stars 
with few photospheric features, and so the importance of 
H I stellar photospheric absorption is diflicult to assess 
in these stars (see Section |473| . 

6.2. Accretion in Transitional Disks 

6.2.1. Equivalent Widths and the Identification of 
Transitional Disks 

In Figure [TT] we show Pf/3 equivalent width (EW) vs. 
Pf /3 luminosity. There is no signiflcant correlation be- 
tween these two variables and thus Pf/3 EW is not a good 
predictor of accretion luminosity. However, we note that 
several transitional disks have notably high Pf/3 EW's 
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Figure 11. Pf/3 equivalent width (EW) plotted against Pf/3 line 
luminosity. Note that EW is not a good predietor for Pf/3 lumi- 
nosity (or therefore accretion luminosity). The dashed line marks 
an EW of —15 A; several transitional disks have extreme values of 
Pf/3 EW. 

(between —20 and —30 A). Since the Pf/3 luminosities of 
these disks are not higher than average, the high EW's 
instead reflect the fact that these disks have reduced con- 
tinuum flux levels in the near-IR. This suggests that Pf/3 
may be a good tracer for some accreting transitional 
disks, and a means to detect these disks with a single 
spectrum — one simple, robust measurement that re- 
quires no knowledge of the stellar parameters or absolute 
flux level. We arbitrarily label all disks with Pf/3 EW 
< — 15 and suggest that these targets may have deple- 
tions of dust in their inner regions. Aside from the tran- 
sitional disks, which are known to have inner disk dust 
depletions, we also note that T Y CrA is a tertiary system 
with an eclipsing binary (e.g.. lVaj|200lD and that Haro 
1-1 is an anomalously fast rotator (jShevchenko &: Herbstl 
[1991) . To our knowledge, no unique properties have been 
discussed for the other disks with high Pf/3 EW. 

It is interesting to note that Pf/3 EW's are not high 
for all transitional disks. In our sample, LkHa 330, HD 
135344 B, IRS 48 and DoAr 44 have typical or even 
slightly low Pf^ EW's. LkHa 330 and HD 135344 B ap- 
pear to have relatively low g as /dust ratios as c ompared 
to other transitional disks (jSalvk et al.l l2009f ) . DoAr 
44 might better be considered a pre-transitional disk 
([Andrews et al.|[2011[ ) — a disk with an inner clearing but 
relatively high near-IR flux, that may be consistent with 
a sm all gap rather than a large clearing ([Espaillat et al.l 
[2001 . Ai id IRS 48 has a ~ 30 AU hole in its gas dis- 
tribution ([Brown et al.ll2012[ ). Therefore, high Pf^ EW 
may be an indicator of only those transitional disks that 
have significantly cleared their inner regions of small dust 
grains but not gas. Thus, a measurement of Pf/3 EW 
might shed light on the process causing the inner disk 
depletion for any given disk, as different clearing scenar- 
ios predict different gas/dust ratios. 

6.2.2. Transitional Disks in Comparison to Other Disks 

Figure [T^ shows the log of Pf/3 luminosity as a func- 
tion of the log of stellar mass. Accretion luminosity and 
accretion rate are known to scale with stellar mass, and 
these results are no exception. Plotting Pf/3 luminosity 
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Figure 12. The top panel shows the log of Pf/3 line luminosity 
plotted against the log of stellar mass and best-fit correlation. The 
lower panel shows the fit residuals. Triangles and diamonds are 
transitional and disk wind targets, respectively; squares are other 
targets. Gray squares represent highly-inclined disks. 



as a function of stellar mass allows us to investigate po- 
tential outliers with high or low accretion rates. 

Transitional disks have on average lower Pf/3 luminosi- 
ties, at a marginally statistically-significant level, with a 
mean log residual Pf/3 luminosity of —0.7 ± 0.3. A low 
Pf/3 luminosity could represent a low accretion luminos- 
ity, or could result if the transitional disks have normal 
accretion luminosities but an anomalously low Pf/3 lu- 
minosity/accretion luminosity ratio. However, we have 
investigated the Pf/3 luminosity /accretion luminosity ra- 
tio for these targets and find that 4 of 5 transitional 
disks with literature accretion luminosities (DM Tau, 
GM Aur, HD 141569 A and TW Hya) have Pf/3 luminosi- 
ties even higher than would be predicted by Equation [T] 
(one, HD 135344 B, has a slightly but not significantly 
lower Pf/3 lumin osity). Therefore, t his st udy is in closer 
agreem ent with iNaiita et al.l (|2007D and lEspaillat et al.l 
(|2012f ) who find slightly lower than average accretion 
rates for transitional disks, than with contrasting stud- 
ies showing no difference between the two types of disks 
([Fang et al. 2009; Sicilia-Aguilar et al. 2010). As sug- 
gested by (Si cilia-Aguilar et al.l I201CD . the contrasting 
study results may be due to the different physical na- 
ture of transitional disks in different regions. However, 
it should be cautioned that the high Pf/3 EW's for tran- 
sitional disks could produce an observational bias in our 
work, making it easier to detect a low accretion rate 
around a transitional, rather than classical, disk. Thus, 
our sample may be missing a number of classical disks 
with low accretion rates. 

Since our sample derives from a large number of clus- 
ters of different ages, it is also worth asking whether we 
might measure spuriously low Pf/3 luminosities for tran- 



sitional disks because they are derived from older clus- 
ters, and there is evidence that accretio n rates decrease 
with stellar age ([Hartmann et al.|ll998D . However, we 
find that the transitional disks are derived from a large 
number of the clusters sampled in this work, and that re- 
gardless of the cluster under consideration, transitional 
disks lie at the low end of the range of measured Pf/3 lu- 
minosities. In addition, the transitional disks are derived 
primarily from clusters with ages < 3 Myr (with the ex- 
ception of TW Hya) , and so are not biased towards older 
ages. 

6.3. Accretion in Disk Wind Sources 

Disk wind targets are a newly identified subset of 
disks defined by their single-peaked near-IR emission 
line flux and spectro-astrometric profiles, which may be 
evidence for the pr esence of a slow-v elocity disk wind 
(|Pontoppidan et al.l[20TTbL iBast et all l201ll ). By com- 
paring average accre tion lumino s ities f or disk wind and 
normal disk targets, IBas t et all (|201lD noted that disk 
wind targets may have higher than average accretion lu- 
minosities. However, Figure [12] shows that the Pf/3 lumi- 
nosities for these targets are actually typical for sources 
in our sample, implying no higher accretion luminosities. 
These targets also have typical ratios of Pf/3 luminos- 
ity to derived accretion luminosity. And, the disk wind 
sources derive from a large number of clusters, implying 
no obvious age-related biases. 

How can these two results be reconciled? We believe 
there are a few possible explana tions. Firstly, the disk 
wind sources in IBast et aD (|2011| ) have somewhat higher 
stellar masses than their comparative sample of "nor- 
mal" disks. Secondly, th e disks in the comparative sam- 
ple in IBast et all ()201lD have lower accretion luminosi- 
ties at a given mass than disks in our sample. Since our 
full sample has not been "vetted" for disk wind targets, 
it is possible that our full sample includes unidentified 
disk wind targets that increase the average, or our sam- 
ple may simply be biased towards higher accr etion rates. 
Unfor tunately, we cannot use the criteria in IBast et al.l 
(|20Tl[ ) to classify all of our targets because of the rela- 
tively lower spectral resolution of NIRSPEC as compared 
to CRIRES. 

The accretion rates for the disk wind targets, derived 
from this work, range from 6xl0~^ to 2x lO~^M0yr~^, 
with most accretion rates in the range of 10~* — 
lO~^Af0yr~^. Thus, the disk wind sources are likely 
capable of supporting mass-loss rates of 10^^° — 
10~^M(7)y r~^, and in some cases 10 ~^M(^yr^^, as sug- 
gested bv iPontoppidan et al.l (|2011bf ). 

Although we find typical accretion rates for disk wind 
sources, it should be noted that one of the primary char- 
acteristics of disk wind targets is t heir high veihng a t 
optical through IR wavelengths (e.g.. lGahm et al.ll2008| ). 
The seemingly contradictory observation of high veiling 
but nor mal accretion rates may be reconciled by the re- 
sults of iGahm et al.l ([2008) , who find that the optical 
veiling in such targets is not a good measure of disk ac- 
cretion, as the derived optical veiling is affected by infill- 
ing of photospheric lines from line emission. 

6.4. Accretion Rates for Embedded Disks 
6.4.1. Classification of Evolutionary Stage 
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The process of star and planet formation is often de- 
constructed into four discrete evolutionary stages. The 
targets in our sample represent primarily, if not en- 
tirely, sources in stages I and II, where stage I sources 
are embedded protostars with still-collapsing envelopes, 
while stage II sources are revealed stars with circumstel- 
lar disks and little or no remaining envelope. For this 
work, the evolutionary class was determined primarily 
through examination of the M-band spectra themselves, 
with stage I disks defined as those showing CO ice ab- 
sorption, and stage II disks those without. This offers the 
significant advantage that our entire sample can be easily 
classified using our own dataset, with a single consistent 
criterion. 

However, it is a known issue that edge-on stage II disks 
can masquerade as evolutionarily-younger systems, es- 
pecially with the use of spectral slope-based classifica- 
tio n schemes (see, fo r example, the detailed discussion 
in lEvans et al.ll2Q09[ ). Our classification scheme is not 
immune to this issue, as CO ices in edge-on disks can in 
principle produce absorption features. However, in prac- 
tice, the production of strong ice absorption features re- 
quires a very specific geometry, with an inclination close 
to the opening angle of the flaring disk. Additionally, the 
upper layers of the disk that intercept the starlight are 
ofte n too warm to contain s ignificant quantities of CO 
ice (jPontoppidan et all I2005D . Thus, it is perhaps not 
surprising that our classification scheme correctly iden- 
tifies some high-inclination disks, including T Tau S and 
VV Ser, as stage II rather than stage I sources. 

Nevertheless, we have made an effort to confirm 
stage I classifications whenever possible using ad- 
ditional diagnostics in the literature, including the 
use of spectral indices, bolometric temperatures, 
the presence of infall or outflow signatures in spec- 
tral lines, and the presence of exten ded envelopes 
as se e n in millimeter- wave n iaps (e. g. iBoogert et aT 



2002; 



2005 



White fc Hillenbrand] 
120091; 



Evans et al.' 
Jorgensen et al. 2009 



I2004t iDoppmann et al 
'van Kempe n et al.l 120091: 

Hcrczcg et al. 20 ll|70ur classifi- 
cations are listed in Table [3l along with accompanying 
notes for sources with ambiguous properties. 

6.4.2. Accretion Rates for Stage I Disks, and Comparison 
with Stage II Disks 

Since stars are believed to accumulate most of their 
mass during their embedded stages, accretion rates for 
stage I disks were initially predicted to be significantly 
larger than for stage II disks. Smooth collapse models 
(jShul II977I ) predicted stage I accretion rates similar to 
simple estimates — for example, accretion rates of or- 
der > lO~^M0yr~-'^ are required to build a solar mass 
star over a ~ 10^ year timescale. Accretion rates in 
stage II disks around solar-mass stars are typically near 

lO~^M0yr~^ (see Tableland references therein), and 
thus stage I accretion rates would be predicted to be at 
least two orders of magnitude higher. 

While there is at least some eviden ce for greater ac- 
cretio n activity in stage I sources (e.g. IDoppmann et al.l 
I2005D . measured accretion rates in stage I disks around 
low-mass stars are found to be near 10~^ Mp ,yr~^ (e.g. 
I White fc Hillenbrandl[200llWhite et alllMl — two or- 
ders of magnitude lower than predicted. This discrep- 
ancy is related to the well-known "luminosity problem" 
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Figure 13. The cumulative distribution function of observed Pf/3 
accretion luminosities for stage I and stage II disks. 



in star formation (jKenvon et al.l I1990D , wherein mea- 
sured bolometric luminosities (which derive primarily 
from the energy of gravitational contraction, and are 
therefore related to mass accretion) are lower than pre- 
dicted from ste ady growth models. Instead, evidence 
is growing (e.g. 'Du nham etahl 12008": 'En och et al] 120091: 
lEvans etah, 2009) that accretion may need to be episodic 
— primarily quiescent, with brief phases of intense ac- 
cretion. 

While direct measurements of accretion rates in stage I 
disks appear to support this hypothesis, such studies in- 
evitably suffer from a bias towards more revealed disks, 
as the most embedded protostars are heavily extincted, 
and cannot be easily studied with either optical or near- 
IR accretion tracers. Although it is difficult to quan- 
tify the many sample-selection biases here, extinction re- 
mains small in the M band even for high Ay, and our 
sample includes several heavily embedded protostars. As 
a point of comparison ,, the s ample of young protostars in 
iWhite fc Hillenbrandl ()2004[ ) all have Ay < 20, while 10 
of 20 stage I disks in our study have Ay > 20. Thus, our 
study provides insight about accretion rates for the most 
embedded, and perhaps evolutionarily youngest, disks. 

Figure [13] shows the empirical cumulative distributions 
of calculated Pf/3 luminosities for stage I and stage II 
disks in our sample. There are small differences between 
the distributions of Pf/3 line luminosities for the two sam- 
ples, with slightly higher Pf/? luminosities in Stage II 
disks. Stage I disks have a mean log(Lpf/3) of —3.7 and 
a standard deviation of 0.8, while stage II disks have a 
mean log(Lpf^) — 3.8 with a standard deviation of 0.9. 
However, a two-sample Kolmogorov-Smirnov test pro- 
duces an associated probability of 30%. Therefore, the 
hypothesis that the stage I and stage II disks draw from 
the same distribution of accretion luminosities cannot be 
rejected, i.e., differences between the two distributions 
are not statistically significant. 

Unfortunately, accretion rate distributions cannot be 
compared directly, as the majority of the embedded disks 
in our sample do not have measured stellar masses or 
radii due to the relative difficulty of obtaining optical 
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spectra for highly embedded stars. However, if average 
stellar masses and radii are similar for the two samples 
our results suggest that accretion rates are quite simi- 
lar for stage I and stage II disks. Spectroscopic analy- 
ses of stage I protostars have suggested that stage I and 
stage II stars have similar average masses, while stellar 
radii are no more than a factor o f ~ 2 larger for stage I 
stars () White fc Hillenbrandll2004 ). For the small number 
(four) of stage I disks with measured stellar masses in our 
sample, the average stellar mass is 1.7Mq, as compared 
to an average of 1.8Mq for our stage II sample. Assum- 
ing at most a factor of two difference in radius between 
stage I and II stars, and assuming no signficantly differ- 
ent selection biases for our two subsamples, our results 
suggest that stage I accretion rates are at most factors of 
a few higher, on average, than they are for stage II disks. 

We can directly measure the accretion rates of four 
stage I targets in our sample. These accretion rates range 
from 1 X 10~^ to 4 X IQ-'^Mgyr-^ Thus, these targets 
have an accretion rate spread of more than an order of 
magnitude, spanning from a value typical for low-mass 
stars to something close to the lO~^M0yr~^ required 
to build a solar mass star in 1 Myr. If we divide the 
currently-measured stellar mass by M, we can obtain a 
lower limit to the time required to form the star assuming 
a steady state accretion rate. We find timescales of 0.7, 6, 
54 and 67 Myr for the four stage I sources — longer than 
estimated stage I lifetimes (Evans et al. 2009). There- 
fore, our data suggest that these sources must at times 
have had significantly higher accretion rates than cur- 
rently observed. 

Making some simple assumptions, we can also investi- 
gate accretion rates for the full stage I sample. For exam- 
ple, assuming M^, = IMq and -R^, = 2Rq for the 16 stage 
I stars without measured stellar parameters, we find a 
mean and standard deviation of log(A/) = —7.1 ± 0.7. 
This value is at least an order of magnitude lower than 
the predictions from steady accretion models. Yet, with 
these assumptions, our sample does include two sources 
with accretion rates of ~ 2 x lO~^M0yr~^. If our sam- 
ple is representative of all stage I sources and homoge- 
neous except for random changes in accretion rate due 
to episodic accretion, and taking into account a detec- 
tion fraction of 60% (see Section [3|), this would imply 
that the stage I sources spend ^ 6% of their lifetime 
with M > lO-^Moyr-^ The large number (40%) of 
Pf/3 non-detections is further evidence for a wide range 
of accretion rates in stage I sources, again consistent with 
an episodic accretion scenario. 

It is also interesting to note that no FU Orionis stars 
in our sample show detectable Pf/3 emission lines. The 
explanation for this is not known, but may have to do 
with the restructuring of the inner disk that occurs dur- 
ing high accretion-rate events fe.g. iZhu et al.|[2009l) . In 
any case, while these sources ma y represent protostars 
with the highest rates of accretion (jHartmann &: KenvonI 
Il996( ). they are necessarily excluded from this analysis 
due to their lack of Pf/3 emission, and thereby bias the 
remaining sample towards lower accretion rates. 

7. CONCLUSIONS 

In this work, we have introduced the use of H i Pf/3 as a 
mass accretion tracer for young stars with protoplanetary 



disks. H I Pf/3 offers several advantages over other accre- 
tion tracers, including being readily observed in heavily- 
extincted disks, requiring minimal correction for photo- 
spheric absorption and being commonly observed along 
with CO fundamental transitions. Using existing mea- 
surements of accretion from the literature, we derive a 
relationship between Pf/3 line luminosity and accretion 
luminosity, and show that this relationship can repro- 
duce measured accretion rates with an accuracy similar 
to that of other commonly-used tracers. Examining our 
large sample of accretion rates, we are further able to 
show that "disk wind" sources appear to have normal 
accretion rates, while our sample of transitional disks 
have slightly lower than average accretion rates. 

We also examine the accretion rates of the stage I 
disks in our sample — a sample that includes signifi- 
cantly more embedded targets than in previous studies 
of this type. We find that stage I and stage II disks have 
statistically indistinguishable Pf/3 luminosities, implying 
similar accretion rates, and that the accretion rates of 
stage I disks are too low to build a stellar mass with qui- 
escent accretion. Our results instead are consistent with 
both observational and theoretical evidence that stage I 
objects experience episodic, rather than quiescent, accre- 
tion. 

The calibration of Pf/3 presented here has allowed for a 
coherent comparison of accretion luminosity in 120 stars. 
While this work has focused exclusively on spectra from 
the NIRSPEC and CRIRES surveys with which we are 
associated, the correlation we derive can be extended to 
a much larger number of protoplanetary disk spectra, in- 
cluding some which may already reside in the NIRSPEC 
or CRIRES archive. We have also focused exclusively 
on disks with clearly detectable Pf/3 emission. However, 
with more extensive observations and focused analysis, 
our correlation can be extended to tenuous disks with 
low accretion rates — disks of utmost importance to our 
understanding of disk dissipation. Finally, as Pf/3 is ob- 
served contemporaneously with CO fundamental emis- 
sion lines, the observation of Pf/3 can be used to ac- 
count for the effects of accretion variability in searches 
for planet-induced CO line variability. We therefore ex- 
pect the utility of Pf/3 as an accretion tracer to extend 
well beyond the work presented here. 
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APPENDIX 

THE SPECIAL CASE OF DOAR 21 

DoAr 21 is known to be an unusual disk, with little 
or no reported excess emission below 7 /im and weak or 
non-existent accretion, but asymmetric extended emis- 
sion beyond ~ 100 AU, and line emission fr om Polycyclic 
Arom atic Hydrocarbons (PAH's) and H2 (jJensen et alj 
120091 and references therein). Observations of DoAr 21 
with the VLB A have revealed t hat DoAr 21 is a 5 mas 
(0.6 AU) binary ()Loinard et all 12008). Detailed studies 
of evolved sources like DoAr 21 can provide insight into 
the process of disk dissipation, as it may represent a case 
of a disk at the latest stages of dissipation. 

DoAr 21 is the only disk in our sample that shows 
both low but detectable Pf/3 equivalent width, and low 
veiling, meaning that the underlying photosphere must 
be accounted for in the calculation of Pf/3 line flux. 
In Figure [TH we show the CRIRES spectrum of DoAr 
21 (Kl) along with a NIRSPEC spectrum of a photo- 
spheric template star — CoKu Tau/4 (K3) — rotation- 
ally broadened and veiled to match the observed spec- 
trum of DoAr 21 longward of 4.68 ^m. Although DoAr 
21 and CoKu Tau/4 do not have the same spectral sub- 
type, we find only small differences in spectra across a 
range of observed spectral types (G2-M0), and both are 
young stars with similarly lowered surface gravities com- 
pared to main-sequence stars. 

The bottom of Figure [14] shows the difference between 
the observed and template spectra, which reveals strong 
and double-peaked Pf/3 and Hue emission lines. As an 
insert in Figure [H} we compare Pf/3 to an observed spec- 
trum of Ha from lJensen et all (|2009f) . The Ha line shown 
here is in fact the difference between Ha lines observed 
on two different nights, and therefore represents the vari- 
able, perhaps flare-related, component to Ha. We find 
that the Pf/3 emission line is broader and more double- 
peaked than the Ha line, suggesting that the Pf/3 may 
not have a flare origin, but instead have an accretion ori- 
gin. (It should be noted, however, that the Hue line flux 
appears unusually high compared to other disks in our 
sample.) 

We also note that the best fit to the observed DoAr 21 
spectrum requires a veiling of 2.6, and that the spec- 
trum is incon sistent with zero veiling, consistent with 
the results of Salv k et "all (|2009f ). Assuming the contin- 
uum flux is emitted entirely by grains at 1500 K with 
an opacity of 500 cm^g~^ (appropriate for sub-micron 
sized grains), the mass in small grains is '^10^^ kg, or 
~ lO""* M J). This is almost certainly a lower limit, as 
there is likely a range of grain sizes and a range of grain 
temperatures, with 1500 K being the approximate maxi- 
mum temperature for solid silicates. An alternative way 
to explain finite veiling would be for the stellar disk to 
contain hotspots with featureless emission spectra. How- 
ever, even in an extreme case of lO'' K hotspots, 54% of 
the stellar disk would need to be covered by hotspots to 
explain the observed veiling. In contrast, results from 
photometric variability studies (jBouvier et al.lflQQSl ) find 
hotspot temperatures of no more than a few hundred to 
a few thousand K hotter than the stellar temperatures, 
and 0.5 — 40% disk covering fractions. Thus, DoAr 21 
may retain both a gas and dust disk at small radii. 
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Figure 14. Top: Spectrum of DoAr 21 (Kl) observed with CRIRES (gray) and a rotationally-broadened, veiled spectrum of CoKu Tau/4 
(K3) observed with NIRSPEC (blaclc). Dotted Unes mark the locations of Pf/3 and Hue. Bottom: Difference spectrum. Inset: Difference 
spectrum centered at 4.6538 /im (gray) compared with observations of Ha from lJensen et aP 1)20091 ) . 
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Table 1 

Log of Observations 



Star 


Date 


Exp 


Vdopa 


Std 




[yyyy mmm dd] 


[«] 


[km s-l] 




NIRSPEC 


AA 'i'au 


2003 Nov 03 


720 


-15 


HR 1620 




2004 Dec 27 


780 


13 


HR 1620 


AB Aur 


2001 Jan 30 


360 


24 


HR 1620 




2001 Aug 07 


180 


-25 


HR 1620 




2002 Jan 03 


240 


13 


HR 1620 


AS 205 N 


2002 Apr 21 


240 


-17 


HR 6175 




2002 Jul 22 


480 


24 


HR 5984 


AS 205 S 


2002 Apr 21 


360 


-17 


HR 6175 


BF Ori 


2009 Dec 28 


240 


6 


HR 1620 


CI Tau 


2008 Oct 10 


480 


-24 


HR 1251 


CK 4 


2006 Jul 07 


600 


3 


HR 7141 


DF Tau 


2005 Dec 21 


240 


10 


HR 1620 


DG Tau 


2001 Oct 26 


225 


-17 


HR 1177 




2005 Dec 12 


120 


6 


HR 1620 




2005 Dec 19 


720 


9 


HR 1620 


DL Tau 


2003 Nov 03 


960 


-15 


HR 1620 


DM Tau 


2003 Nov 03 


720 


-14 


HR 1620 




2006 Oct 09 


1380 


-24 


HR 1251 


DoAr 44 


2002 Jul 22 


480 


22 


HR 5984 




2004 Jul 23 


720 


23 


HR 6175 




2005 Apr 26 


480 


-17 


HR 6070 


DO Tau 


2005 Dec 18 


240 


8 


HR 1620 




2005 Dec 20 


480 


9 


HR 1620 


DP Tau 


2006 Dec 29 


480 


13 


HR 1620 


DR Tau 


2002 Dec 17 


480 


7 


HR 1620 




2005 Dec 18 


240 


7 


HR 1620 


DS Tau 


2007 Dec 26 


420 


10 


HR 1620 


FS CMa 


2002 Dec 17 


240 


-6 


HR 1620 


FT Tau 


2008 Dec 10 


900 


5 


HR 1620 


FZ Tau 


2008 Dec 08 


480 


3 


HR 1620 


GM Aur 


2004 Dec 27 


780 


10 


HR 1620 


GQ Lup 


2006 Apr 15 


840 


-18 


HR 5812 


GV Tau S 


2008 Dec 10 


630 


5 


HR 1620 


GY 224 


2006 Jul 07 


420 


17 


HR 7141 


Haro 1-1 


2005 Apr 25 


840 


-17 


HR 6175 




2009 Jul 13 


240 


19 


HR 5993 




2009 Jul 13 


480 


20 


HR 5993 


Haro 1-4 


2009 Jul 13 


300 


20 


HR 5993 


Haro 6-13 


2008 Dec 08 


480 


3 


HR 1620 


HD 135344 B 


2005 Apr 24 


480 


-11 


HR 5712 




2006 Jul 06 


240 


20 


HR 5812 


HD 141569 A 


2002 Apr 21 


210 


-13 


HR 6175 




2008 Jul 09 


480 


22 


HR 5685 


HD 142666 


2002 Jul 21 


240 


24 


HR 5984 


HD 143006 


2005 Apr 26 


480 


-14 


HR 6175 


HD 144432 S 


2002 Jul 21 


480 


23 


HR 5984 




2005 Apr 26 


480 


-15 


HR 6175 


HD 150193 


2002 Jul 22 


480 


22 


HR 5984 


HD 163296 


2001 Aug 08 


150 


21 


HR 6581 




2002 Apr 21 


240 


-25 


HR 6378 




2002 Jul 22 


480 


15 


HR 7141 


HD 169142 


2002 Jul 21 


420 


11 


HR 6175 




2006 Jul 07 


240 


5 


HR 7141 


HD 179218 


2006 Oct 10 


180 


23 


HR 7235 


HD 190073 


2002 Jul 22 


480 


-2 


HR 7141 




2005 Sep 22 


600 


22 


HR 7724 


HD 244604 


2009 Dec 28 


240 


7 


HR 1620 


HD 35187 


2009 Dec 28 


420 


8 


HR 2714 


HD 37258 


2009 Dec 28 


240 


6 


HR 1620 


HD 37357 


2009 Dec 28 


240 


6 


HR 1620 


HD 37411 


2009 Dec 28 


420 


6 


HR 2714 


HD 37806 


2009 Dec 28 


240 


5 


HR 2714 


HD 38120 


2009 Dec 28 


240 


5 


HR 2714 


HD 50138 


2009 Dec 28 


240 


-4 


HR 2714 


HD 58647 


2008 Dec 08 


240 


-16 


HR 1620 




2009 Dec 28 


480 


-8 


HR 2714 


HL Tau 


2001 Oct 25 


200 


-18 


HR 1380 




2002 Dec 17 


240 


8 


HR 1620 


IM Lup 


2006 Jul 06 


420 


18 


HR 5812 


IRAS 03259-1-3105 


2002 Dec 18 


240 


14 


HR 1177 


IRAS 03301-1-3111 


2001 Oct 26 


200 


-12 


HR 1177 




2002 Dec 16 


720 


13 


HR 1177 




2002 Dec 17 


360 


13 


HR 1177 
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Table 1 — Continued 



Star 


Date 


Exp 


Vdop^ 


Std 




lyyyy nimm aaj 


L=J 


^km s j 




IRAS 03445-1-3242 


2002 Dec 16 


480 


11 


HR 1177 




2001 Oct 26 


300 


-14 


HR 1177 


IRAS 042394-2436 


2005 Dec 12 


360 


6 


HR 1620 


IRS 51 


2002 Apr 21 


300 


-19 


HR 6175 


IRS 60 S 


2002 Jul 21 


480 


22 


HR 6175 


LkHa 25 


2006 Apr 16 


480 


28 


HR 2714 


LkHa 198 


2002 Jul 22 


480 


-18 


HR 7945 


LkHa 220 


2001 Jan 30 


600 


9 


HR 1620 


LkHa 271 


2009 Dec 28 


780 


18 


HR 1620 


LkHa 326 


2002 Jul 22 


960 


-25 


HR 1177 




2002 Dec 16 


600 


13 


HR 1177 


LkHa 327 NE 


2002 Dec 16 


240 


13 


HR 1177 




2002 Dec 16 


480 


13 


HR 1177 


LkHa 327 SW 


2002 Dec 16 


240 


13 


HR 1177 




2002 Dec 16 


480 


13 


HR 1177 


LkHa 330 


2002 Dec 16 


480 


11 


HR 1177 




2003 Nov 02 


360 


-11 


HR 1620 




2003 Nov 03 


960 


-10 


HR 1177 




2004 Dec 27 


480 


17 


HR 1177 


LkHa 348 


2002 Apr 21 


120 


-26 


HR 7141 




2009 Jul 12 


240 


5 


HR 7141 




2002 Apr 21 


300 


27 


HR 2763 




2002 Dec 16 


240 


-7 


HR 2034 




2004 Dec 27 


480 


-1 


HR 2763 


MWC 480 


2001 Jan 30 


480 


24 


HR 1620 




2001 Aug 08 


360 


-25 


HR 1620 




2002 Dec 17 


480 


4 


HR 1620 


MWC 758 


2002 Apr 21 


120 


24 


HR 2763 




2002 Dec 16 


240 





HR 2034 


MWC 1080 


2002 Jul 22 


480 


-16 


HR 7945 


RNO 19 


2005 Dec 19 


900 


13 


HR 1177 


RNO 90 


2002 Jul 21 


600 


22 


HR 5984 


RNO 91 


2002 Jul 21 


600 


22 


HR 5984 




2003 Jul 10 


180 


18 


HR 6175 


ROX 43A 


2002 Apr 21 


300 


-19 


HR 6175 


RR Tau 


2005 Dec 18 


180 





HR 1620 


RXJ1842.9-3532 


2005 Apr 26 


480 


-26 


HR 7121 


RY Tau 


2002 Dec 17 


720 


9 


HR 1620 


SR 4 


2005 Apr 25 


960 


-17 


HR 6175 


SU Aur 


2002 Dec 16 


240 


4 


HR 1620 




2003 Nov 02 


240 


-18 


HR 1620 




2003 Nov 03 


960 


-17 


HR 1620 


T Tau S 


2002 Jan 03 


360 


18 


HR 1620 




2002 Dec 16 


240 


9 


HR 1177 




2002 Dec 17 


240 


9 


HR 1177 




2004 Dec 27 


480 


15 


HR 1620 


T Tau N 


2002 Jan 03 


360 


18 


HR 1620 




2002 Dec 16 


240 


9 


HR 1177 




2002 Dec 17 


240 


9 


HR 1177 




2004 Dec 27 


480 


15 


HR 1620 


TW Hya 


2001 Jan 29 


1020 


-19 


HR 4494 




2001 Jan 30 


1170 


-19 


HR 4494 




2002 Apr 21 


570 


11 


HR 4313 




2002 Dec 17 


1020 


-24 


HR 1620 




2004 Dec 27 


300 


-24 


HR 4313 




2005 Apr 25 


960 


12 


HR 4494 




2005 Dec 18 


960 


-24 


HR 4494 


UY Aur A 


2006 Dec 29 


420 


11 


HR 1620 




2007 Oct 30 


420 


-18 


HR 1620 


UY Aur B 


2007 Oct 30 


180 


-18 


HR 1620 


UX Ori 


2005 Dec 12 


300 


3 


HR 1620 




2005 Dec 19 


720 


6 


HR 1620 


VSSG 1 


2008 Jul 10 


120 


19 


HR 5993 


VV Ser 


2001 Aug 06 


600 


16 


HR 7141 




2003 Jul 09 


1380 


4 


HR 7141 


V1685 Cyg 


2005 Sep 21 


720 


9 


HR 7546 


V1977 Cyg 


2005 Sep 21 


600 


6 


HR 7546 


Wa Oph 6 


2005 Apr 26 


240 


-18 


HR 6175 


WL 16 


2002 Apr 21 


240 


-19 


HR 6378 


WL 22 


2005 Apr 25 


330 


-17 


HR 6175 


WW Vul 


2002 Jul 22 


480 


1 


HR 7141 




2005 Sep 22 


660 


19 


HR 7724 


CRIRES 


AA Tau 


2007 Oct 13 


600 


-23 


HR 1178 


AS 205 N 


2007 Apr 22 


1200 


-16 


HR 4757 
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Table 1 — Continued Table 2 — Continued 



Earth-induced Doppler shift towards the source position. 



Table 2 

Photometry 



Star 



M 



4.5=^ 



5.8=^ Refs 



AA Tau 
AB Aur 
AS 205 N 
AS 205 S 
AS 209 
BF Ori 
CI Tau 
CK 4 

CrA IRS 2 
CV Cha 
CW Tau 
DF Tau 
DG Tau 
DL Tau 
DM Tau 
DoAr 21 
DoAr 24 E S 
DoAr 24 E N 



6.86 
2.70 
4.18 
5.31 



6.77 6.35 



6.27 
7.49 



5.05 

5.94 
7.56 



12 


151 


16 




17 




17 




w 



24l 

El 



4.03 






5.09 






5.17 




5.28 


4.63 




4.19 


6.04 


6.12 


5.66 




9.19 


9.25 


5.51 


5.87 


5.51 



8lT2 



IRS 44 E 

IRS 48 
IRS 51 
IRS 60 S 
IRS 63 
LkHa 25 
LkHa 198 



21<= 



5.58 
6.91 



5.12 
5.73 
6.10 



otar 


Date 


Exp 


Vclop 


ota 


otar 


1\/T 
iVl 


4.0 


o.o 




lyyyy mmm dd] 


r 1 
[s] 


n —11 
[km s '-J 




DoAr 44 




6.30 


5.94 


AS 209 


2008 May 02 


480 


— 15 


HR 6084 


DO Tau 


5.27 


5.65 


5.20 


CrA TRS 2 


2007 Apr 26 


480 


—26 


HR 7236 


DP Tau 




6.91 


6.37 


CV Cha 


2009 Jan 04 


1200 


—9 


HR 4467 


DR Tau 


4.80 




4.50 


CW Tail 


2008 Dec 31 


720 


16 


HR 2451 


DS Tau 




6.91 


6.62 


DF Tau 


2009 Jan 02 


720 


16 


HR 2550 


EC 90 N 


5.76 






DG Tau 


2007 Oct 16 


720 


—21 


HR 2550 


EC 90 S 


4.50 






DoAr 21 


2008 May 01 


720 


— 14 


HR 7236 


Ehas 29 


2.40 






DoAr 24 E S 


2007 Sep 03 


600 


29 


HR 7913 


EX Liin 


3 ggb 






DoAr 24 E N 


2007 Sen 03 


600 


29 


HR 791 3 


FS Cm A 








DoAr 44 


2007 Sen 05 


960 


29 


HR 6084 


FT Tau 




7.15 


6.88 


DR Tau 


2007 Oct 1 1 


480 


—24 


HR 31 1 7 


FZ Tau 




5.69 


5.21 


EC 90 N 


2007 Apr 24 


600 


—24 


HR 5531 


CM Aiir 




7.88 


7.70 


EC 90 S 


2007 Apr 24 


600 


—24 


HR 5531 


GO Liin 






5.64 


Elias 29 


2008 Aug 09 


360 


27 


HR 7235 


GV Tau S 








EX Lup 


2008 Apr 28 


480 


— 13 


HR 6084 


GY 224 




6.92 


6.28 


GO Liin 


2007 Apr 22 


1200 


— 15 


HR 4757 


Harn 1-1 


8.03 




9.80 


Haro 1-4 


2008 Anr 30 


720 


— 14 


HR 6074 


Harn 1-4 






6.24 




2007 Anr 23 


2400 


— 11 


HR 5531 


Harn 6-1 3 


5.65 






HD 142527 


2008 All? 07 


360 


25 


HR 5984 


HD 1 35344 R 


4.69 






HD 144432 S 


2008 Aiie- 03 


2880 


26 


HR 6508 


HD 1 41 56Q A 


6.58 






HD 144965 


2008 Anr 30 


720 


— 13 


HR 6084 


HD 142527 


3.50 






HH 100 


2007 Sen 03 


600 


24 


HR 6508 


HD 142666 


4.69 






HL Tau 


2007 Oct 12 


600 


—22 


HR 838 


HD 143006 


5.34 






IRS 43 


2008 All? 06 


960 


26 


HR 6175 


HD 144432 S 


4.90 






IRS 44 E 


2008 Aug 07 


960 


26 


HR 6175 


HD 144965 








IRS 48 


2008 May 03 


240 


— 13 


HR 6084 


HD 1 501 93 


3.93 






IRS 51 


2008 Aiie- 03 


1200 


25 


HR 6175 


HD 163296 


3.14 






TRS 63 


2007 Apr 26 


600 


— 17 


HR 5531 


HD 169142 


5.57 






LkHa 330 


2008 Dec 30 


720 


17 


HR 1149 


HD 179218 


4.18 






LLN 8 


2007 Oct 12 


720 


— 12 


HR 3117 


HD 190073 








LLN 17 


2007 Oct 14 


1080 


— 11 


HR 31 1 7 


HD 244604 


5.37 






RNO 91 


2007 Apr 24 


600 


— 17 


HR 5531 


HD 35187 








RU Lup 


2007 Apr 27 


600 


— 13 


HR 5531 


HD 37258 








S CrA A 


2007 Apr 23 


480 


—27 


HR 5531 


HD 37357 








S CrA B 


2007 Apr 23 


480 


—27 


HR 5531 


HD 37411 


5.81 






Sz 68 


2007 Apr 27 


600 


— 12 


HR 5531 


HD 37806 


3.79 






T CrA 


2007 Apr 24 


600 


—27 


HR 5531 


HD 38120 








TW Hya 


2007 Apr 22 


1200 


11 


HR 4757 


HD 50138 








TY CrA 


2007 Apr 25 


600 


-26 


HR 7236 


HD 58647 


3.50 






VSSG 1 


2008 Aug 07 


960 


26 


HR 6175 


HH 100 








VV CrA S 


2007 Apr 25 


600 


-26 


HR 7236 


HL Tau 


4.09 






VV CrA N 


2007 Apr 25 


600 


-26 


HR 7236 


IM Lup 


6.55 


7.28 


6.43 


VW Cha 


2009 Jan 04 


720 


-8 


HR 4467 


IRAS 032594-3105 




6.74 


4.05 


VZ Cha 


2009 Jan 04 


720 


-9 


HR 4467 


IRAS 03301-1-3111 




6.53 


5.65 


Wa Oph 6 


2008 May 02 


480 


-15 


HR 6084 


IRAS 03445-1-3242 




5.86 


4.80 


WL 12 


2007 Sep 02 


960 


29 


HR 5812 


IRAS 04239-1-2436 




6.32 


5.38 


WL 22 


2008 Aug 04 


1200 


26 


HR 6175 


IRS 43 




5.00 


4.01 



Refs 



3.64 
5.17 
5.45 



El 


121151 


8 
8 


12| 


8 

g 




9 




9 




7 





5l26l 



1^ 

M 

16 

18 
16 
16 
IS 

IS 

M 

16 
16 
16 



m 

111221 



El 

I 

23] 

i 

5] 
SI 
SI 
9] 
26] 



LkHa 220 


7.00 






26] 


LkHa 271 




8.22 


8.02 


S] 


LkHa 326 




6.97 


6.49 


S] 


LkHa 327 NE 










LkHa 327 SW 










LkHa 330 




5.77 


5.48 


El" 


LkHa 348 








LLN 8 


4.01 




■■■ 


H 


LLN 17 


4.70 




... 


H 


MCW 442 


3.50 








MWC 480 










MWC 758 


4.47 






El' 


MWC 1080 


2.40 






m 


RNO 19 




6.78 


6.22 


g 


RNO 90 






4.43 


SI 


RNO 91 










ROX 43A 




5.88 


5.29 




RR Tau 


5.50 








RU Lup 

RXJ1842.9-3532 
RY Tau 


5.38 
3.38 


7.30 


5.12 


5|25] 


ft 



20 
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Star 


M 


4 5a 




Refs 


S CrA A 


4.15 






m 


S CrA B 


5.33 






m 


SR 4 


6.90 


6.19 


5.86 


im 


SU Aur 


4.26 




4.64 


8|12| 


Sz 68 




5.25 


4.60 


m 


T CrA 


5.18 








T Tau S 










T Tau N 










TW Hya 




7.07 


6.93 




TY CrA 


5.83 






i 


UY Aur A 










UY Aur B 










UX Ori 


5.32 




■ ■■ m 


VSSG 1 




5.62 


5.20 


m 


VV CrA S 










VV CrA N 










VV Ser 


4.40 


4.63 


3.83 


5|10| 


VW Cha 










VZ Cha 










V1685 Cyg 










V1977 Cyg 










Wa Oph 6 




5.72 


5.79 


p 


WL 12 


5.76 






91 


WL 16 




4.90 


3.41 


5] 


WL 22 




6.42 


5.79 


5] 


WW Vul 









References. — (1) Bouvier et al. 1992 (2) Carpenter et al. 
2008 (3) Cieza et al. 2007 (4) Ducati et al. 2002 (5) Evans et al. 
2003 (6) Evans et al. 2009 (7) Goto et al. 2011 (8) Hartmann et 
al. 2005 (9) Herczeg et al. 2011 (10) Hillenbrand et al. 1992 (11) 
Hughes et al. 1994 (12) Kenyon & Hartmann 1995 (13) Lada et al. 

2006 (14) Liseau et al. 1992 (15) Luhman et al. 2006 (16) Malfait 
et al. 1998 (17) McCabe et al. 2006 (18) Men'n et al. 2004 (19) 
Merm et al. 2008 (20) Padgett et al. 2006 (21) Padgett et al. 2007 
(22) Pinte et al. 2008 (23) RebuU et al. 2010 (24) Robitaille et al. 

2007 (25) Schegerer et al. 2009 (26) Wright et al. 2010 

Spitzer-IRAC filters. 
^ Although other photometric measurements are available for this 
source, we use the spectroscopically-derived measurements from 
I Goto et al.l II2011I ) as EX Lup was in outburst during the observa- 
tions presented here, and thus several magnitudes brighter than 
in its quiescent state. 

Corrected for observed relative fluxes of IRS 44 E and IRS 44 
W from lHerczee et al.l 1I20T1I') . 
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Table 3 

Stellar and disk parameters 



Dtar 


dist 


op 1 


1 ^ 


T 






Ay 


AA Tau 


140 


K7 


4060 


0.71 


1.74 


0.53 


0.74 


AB Aur 


140 


AO 


9840 


48 




2.40 


0.5 


AS 205 N 


160 


K5 


4450 


7.1 




1.50 


2.9 


AS 205 S 


160 


M3 


3450 


2.19 




0.30 


2.1 


AS 209 


160 


K5 


4395 


2.5 


2.40 


1.40 


1.15 


BF Ori 


400 


A5-F6 




27 


2.70 


2.10 


2.4 


CI Tau 


140 


K7 


4060 


0.84 




0.77 


1.8 


CK 4 


260 


K3 


4730 








12.8 


CrA IRS 2 


130 


K2 


4900 








20 


CV Cha 


215 


G8 


5451 


8.0 


3.20 


2.10 


1.67 


CW Tau 


140 


K3 


4730 


0.76 




1.06 


2.21 


DF Tau 


140 


MO. 5 


3470 


1.97 


3.37 


0.27 


0.45 


DG Tau 


140 


K6 


4200 


0.28 


1.00 


0.85 


1.3 


DL Tau 


140 


K7 


4060 


1.16 




0.76 


2.0 


DM Tau 


140 


Ml 


3705 


0.32 




0.62 


1.1 


DoAr 21 


122 


Kl 


5080 


11.7 




1.80 


6.0 


DoAr 24 E S 


120 


KO 


5250 


8.78 




0.47 


5.7 


DoAr 24 E N 


120 












5.7 


DoAr 44 


120 


K3 


4730 


1.4 


1.75 


1.30 


2.2 


DO Tau 


140 


MO 


3850 


1.0 


2.25 


0.37 


2.3 


DP Tau 


140 


MO 


3850 


0.18 




0.57 


0.6 


DR Tau 


140 


K7 


4060 


1.7 




0.40 


1.7 


DS Tau 


140 


K5 


4350 


0.57 


1.36 


0.87 


0.34 


EC 90 N 


415 


M5 


3270 








4.0 


EC 90 S 


415 


Gl 


5900 








30 


Elias 29 


120 












9.8 


EX Lup 


140 


MO 


3802 


0.39 









FS CmA 


500 




21878 


3162 






1.0 


FT Tau 


140 


Cont 


3890 











FZ Tau 


140 


MO 


3918 


0.51 




0.70 


2.72 


CM Aur 


140 


K7 


4060 


0.74 


1.78 


0.52 


0.31 


GQ Lup 


100 


K7 


4060 


0.8 




0.80 


0.5 


GV Tau S 


140 


K4 


4500 


4.9 






30 


GY 224 


120 


M4 


4635 


3.7 




0.91 


34.2 


Haro 1-1 


120 


K5-K7 


4205 








1.67 


Haro 1-4 


120 


K4-K6 


4350 








3.8 


Haro 6-13 


140 


MO 


3800 


2.11 




0.52 


11.9 


HD 135344 B 


84 


F3 


6590 


8.0 




1.50 


0.3 


HD 141569 A 


99 


B9.5 


9550 


22.9 


1.80 


2.20 


0.1 


HD 142527 


200 


F6 


6270 


69 




3.50 


0.7 


HD 142666 


145 


A8 


7590 


17 


2.40 


2.00 


0.8 


HD 143006 


145 


G5 


5884 


2.5 






1.63 


HD 144432 S 


145 


A9 


7410 


14.8 


2.30 


2.00 


0.2 


HD 144965 


190 


B5 













HD 150193 


203 


Al 


8970 


36.1 


2.50 


2.20 


1.5 


HD 163296 


122 


Al 


9450 


36 




2.30 


0.3 


HD 169142 


145 


A5 


8260 


15 




1.80 


0.3 


HD 179218 


240 


B9 


10600 


500 




4.00 


1.8 


HD 190073 


767 


A2 


9500 


471 


8.00 


5.10 


0.5 


HD 244604 


336 


AO 


9549 


97.7 




3.05 


1.69 


HD 35187 


150 


A2 


9120 


27 




2.30 


0.65 


HD 37258 


510 


A2 


8970 








0.3 


HD 37357 


240 


A2 


10495 


29.5 




2.40 


0.35 


HD 37411 


510 


B9 


10500 








0.52 


HD 37806 


470 


B9 


8912 


138 




3.58 





HD 38120 


460 


B9 


10471 


123 




3.12 


0.67 


HD 50138 


500 


B6 


1148 








0.2 


HD 58647 


543 


B9 


10500 


911 


9.10 


6.00 


1.0 


HH 100 


130 


K7 


4060 


3.1 


3.60 


0.40 


30 


HL Tau 


140 


K5 


4350 


0.9 




0.55 


3.2 


IM Lup 


140 


MO 


3918 


1.3 






0.98 


IRAS 03259-f 3105 


250 












5.9 


IRAS 03301-1-3111 


250 














IRAS 03445-1-3242 


350 




5750 


9.40 






30 


IRAS 04239-1-2436 


140 






0.15 






33 


IRS 43 


120 


K5 


4400 








9.8 


IRS 44 E 


120 


K5 


4400 








9.8 


IRS 48 


120 


AO 


9000 


14.3 




2.00 


11.5 


IRS 51 


120 


K6 


5200 


1.1 






37.6 


IRS 60 S 


120 


K2 










9.0 


IRS 63 


120 












9.8 


LkHa 25 


800 


AO 










0.65 


LkHa 198 


600 


A5 










4.5 


LkHa 220 


1050 














LkHa 271 


250 


K3-K5 










0.9 



Rofs 



Stage 



Ref Inst. ^ 



I27I60I 

imT 

[52] 

[Toil 

I30I38I55I 

I5I43I 
I38I104I 



6113140 



9151163 



mm 




I14I70I 
[Toil 




I32I69I 
I48I78I102I 

I32I69I 
I21I50I76I 

I32I69I 



[32] 
I69I73I95I 
Ill(i5l67l 
[32] 
I65I67I106I 
I106I97I 



I8I65I106I 
I65I6VI9VI 
I46I67I99I 



111106198 



l-'4|j-.|.:l| 



[861 
||."| I 



II 
II 
II 
II 
II 
II 
II 
IF 
I 
II 
II 
II 
II 
II 
lit 
II 
II 
II 
lit 
II 
II 
II 
II 
I 
I 
I 
II 
II 
II 
II 
lit 
II 
I 
I 

II 
II 

Ild 
lit 
lit 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 
I 

P 
II 

I 
Ilf 

I 

I 

I 

I 

lit s 
I h 
II 
I 

II 
II 
II 
II 



[44] 
[H 



[27] 
I24I44I 

[g 



[H 



22 
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Star 


dist 


SpT 


T* 










LkHa 326 


250 


G-MO 












LkHa 327 NE 


250 












3.83 


LkHa 327 SW 


250 












3.83 


LkHa 330 


250 


G3 


5800 


16 




2.50 


1.8 


LkHa 348 


260 














LLN 8 


400 






642 




5.00 


20 


LLN 17 


700 












20 


MCW 442 


800 


B6 


14125 


1550 






1.2 


MWC 480 


137 


A5 


8200 


21.9 




1.97 


0.33 


MWC 758 


204 


A3 


8130 


30 




2.00 


0.66 


MWC 1080 


2200 


BO 


30200 


5750 






8.4 


RNO 19 


250 












3.9 


RNO 90 


120 


G5 


5662 


5.7 




1.60 




RNO 91 


120 












9.0 


ROX 43A 


120 


GO 


6030 








3.5 


RR Tau 


2103 


AO 


10000 


781 


9.30 


5.80 


1.1 


RU Lup 


140 


K7 


4060 


0.42 


1.30 


1.00 


0.1 


RXJ1842.9-3532 


140 


K2 


4995 


1.12 






1.03 


RY Tau 


140 


Gl 


5945 


9.6 


2.90 


2.00 


2.2 


S CrA A 


130 


K3 


4800 


2.29 




1.50 


1.0 


S CrA B 


130 


MO 


3800 


0.76 




0.60 


1.0 


SR 4 


120 


K6 


4197 


1.5 




0.52 


1.9 


SU Aur 


140 


Gl 


5945 


7.80 


2.60 


1.70 


0.9 


Sz 68 


159 


K2 


4900 


14.49 




2.50 


4.0 


T CrA 


130 


FO 


7244 


7.90 






2.1 


T Tau S 


140 










2.50 


35 


T Tau N 


140 


KO 


5248 


7.3 




2.11 


1.46 


TW Hya 


56 


K7 


4060 


0.17 


0.83 


0.77 





TY CrA 


130 


B9 










1.75 


UY Aur A 


140 












2.1 


UY Aur B 


140 












2.1 


UX Ori 


450 


A3 


8670 


39 




2.30 


0.3 


VSSG 1 


120 


MO 


3890 


1.5 




0.52 


0.97 


VV CrA S 


130 


Kl 










2.4 


VV CrA N 


130 


Kl 










2.4 


VV Ser 


260 


B9 


10600 


85 




3.10 


3.1 


VW Cha 


150 


K5 




3.37 


2.75 


0.60 


2.39 


VZ Cha 


150 


K7 


4060 


0.46 


1.30 


0.78 


0.47 


V1685 Cyg 


1000 












3.16 


V1977 Cyg 


826 


B8 


11000 


207 


4.00 


3.50 


1.74 


Wa Oph 6 


120 


K 




0.67 






3.5 


WL 12 


120 


K7 


4000 








9.8 


WL 16 


120 


K5 


4467 


4.57 




1.02 




WL 22 


830 












9.8' 


WW Vul 


696 


A3 


8970 


50 


2.90 


2.50 


1.0 



Rcfs 



Stage 



Ref Inst. 



601106 
601106 
I15I60I 





18|68 




C 




62|68 




C 




[891 




N 


I72I73I95I n 


N 



I33I41I87I 

EH 

I32I69I 





m 

78] 




21 

m 




m 


I9I12I75I 



EI] 
I1I5TI106I 



I28I40I 
I28I40I 
I91I106I 
I69I73I106I 
I2I34I100I 
\M 
[75] 

M 
I32I69I 
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Star dist SpT T^, _R* Ay Rcfs Stage Ref Inst. 



References. — (1) Acke et al. 2005 (2) Andrews & Williams 2007 (3) Andrews ct al. 2011 (4) Barsony ct al. 2005 (5) Basri & Batalha 
1990 (6) Beckwith et al. 1990 (7) Blondel et al. 2006 (8) Boersma et al. 2008 (9) Bontemps et al. 2001 (10) Boogert et al. 2002 (11) Borges 
Fernandes et al. 2009 (12) Bouvier & Appenzeller 1992 (13) Briceiio et al. 2002 (14) Brown et al. 1995 (15) Brown et al. 2007 (16) Brown 
et al. 2012 (17) Calvet et al. 2004 (18) Caratti o Garatti et al. 2004 (19) Chini 1981 (20) Cohen & Kuhi 1979 (21) Correia et al. 2006 (22) 
de Lara et al. 1991 (23) Doppmann et al. 2003 (24) Doppmann et al. 2005 (25) Duehene et al. 2004 (26) Dunkin & Crawford 1998 (27) 
Evans et al. 2009 (28) Feigelson et al. 1993 (29) Friedemann et al. 1993 (30) Furlan et al. 2006 (31) Furlan et al. 2009 (32) Garcia-Lopez 
et al. 2006 (33) Chez, Neugebauer & Matthews 1993 (34) Gras- Velazquez & Ray 2005 (35) Greene & Meyer 1995 (36) Gredel et al. 1997 
(37) Guenther et al. 2007 (38) Gullbring et al. 1998 (39) Hartigan & Kenyon 2003 (40) Hartmann et al. 1998 (41) Herbst & Warner 1981 
(42) Herczeg et al. 2004 (43) Herezeg et al. 2008 (44) Herezeg et al. 2011 (45) Hernandez et al. 2004 (46) Hernandez et al. 2005 (47) Heyer 
et al. 1990 (48) Hillenbrand et al. 2008 (49) Hodapp 1994 (50) Hughes et al. 1994 (51) Jensen et al. 2009 (52) Johns-Krull et al. 2000 
(53) Johns-Krull & Gafford 2002 (54) Keller et al. 2008 (55) Kenyon & Hartmann 1995 (56) Kessler-Silacci et al. 2006 (57) K ohler et al. 
2008 (58) Koresko et al. 1997 (59) Lada et al. 2006 (60) Lahuis et al. 2007 (61) Levreault 1985 (62) Liseau et al. 1992 (63) Loinard et 
al. 2008 (64) Luhman & Rieke 1999 (65) Malfait et al. 1998 (66) Mannings & Sargent 2000 (67) Manoj et al. 2006 (68) Massi et al. 2006 
(69) Mendiguti'a et al. 2011 (70) Miroshnichenko 2007 (71) Monnier et al. 2005 (72) Montesinos et al. 2009 (73) Mora et al. 2001 (74) 
Myers et al. 1987 (75) Natta et al. 2006 (76) Neckel et al. 1980 (77) Nisini et al. 2005 (78) Pascucci et al. 2007 (79) Ferryman et al. 1997 
(80) Pontoppidan et al. 2007 (81) Prato et al. 2003 (82) Prato et al. 2009 (83) Riee et al. 2006 (84) Seperuelo Duarte et al. 2008 (85) 
Shevchenko & Herbst 1998 (86) Shevchenko et al. 1999 (87) Simon et al. 1992 (88) Steele et al. 1999 (89) Stelzer et al. 2009 (90) Straizys 
et al. 1996 (91) Strom et al. 1972 (92) Testi et al. 1998 (93) Thompson et al. 1998 (94) Valenti et al. 1993 (95) Valenti et al. 2003 (96) 
van den Ancker et al. 1998 (97) van Boekel et al. 2005 (98) van Leeuwen 2007 (99) Vieira et al. 2003 (100) Walter 1986 (101) Watson et 
al. 2009 (102) White et al. 2007 (103) White &; Hillenbrand 2004 (104) White & Chez 2001 (105) Wichmann et al. 1998 (106) Yudin 2000 

^ lit refers to known transitional disks with inner disk cavities. 

^ Instrument. N=NIRSPEC, C=CRIRES. 

CK 4 (SSTc2d J182958. 19+011521.8) has a weak ice absorption feature in NIRSPEC spectra, but its spectral slope, bolometric 
temperature and lack of extended millimeter emission suggest that it is a stage II object (Evans et al. 2009). 

Haro 6-13 is often classified as a borderline stage I/II object l ITeixeira fc Emerson 1999 : White & H illenbrandll200^ ) . 
° HL Tau is often classified as a borderline stage I/II object lITeixeira Ik Emerson|[T99SI : IWhite fc Hillenbrandll2004l) , but is classified as 
stage I bv lHerczeg et aLl I I2011I ) due to the presence of an extended envelope, cold compact gas and molecular outflow. 

This source is borderline stage I/II based on its "flat" spectral slope and bolometric temperature, and the presence of an extended 
envelope ijEvans et al.|[2009h . However, its M-band spectrum shows no evidence for CO ice, and we adopt a stage II classification. 
^ IRS 48 has a weak CO ice feature in its CRIRES spectrum, but appears to be an evolved stage II disk, perhaps transitioning to stage 
III l IBrown et al.ll20T^ . 

^ IRS 51 is a flat spectrum source with a bolometric temperature in the Class II range, but with evidence of an extended envelope 
l lEvans et al.|[2009t) . We adopt a stage I classification due to its prominent M-band ice feature. 
' Assuming the properties of SSTc2d -I-J162659.10- 243503.3 l lEvans et al.|[2009D . 
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Table 4 

Derived Pf/3 Accretion Parameters; EW's of Other Accretion Tracers 



iDtar 


TAT 

[A] 


1 r a 
log Lpfp 


iOg -Lace 
r r 1 


Inn- l\'/r^ 
log Ivl 

[Moyr ^] 


rj vv J 
r A 1 
[A] 


Hj VV \^rd,p J 
r A 1 

[A] 


xxeis 


AA Tau 


—2.21 


—5.21 


— 1.43 


—8.31 


—0.59 


—2.58 


2 


AB Aur 


—6.75 


—3.09 


0.50 


—6.90 


—5.50 




3 


AS 205 N 


-3.59 


-3.83 


-0.18 


-7.10 








AS 205 S 


—5.29 


—4.10 


—0.42 


—6.68 








AS 209 


—4.41 


—4.02 


—0.35 


—7.52 








BF Ori 


—4.74 


—3.69 


—0.05 


-7.34 


—3.60 






CI Tau 


—8.88 


—4.37 


—0.66 


—7.68 


—4.20 


—7.00 


5 


CK 4 


—12.67 


—4.04 


—0.37 










CrA IRS 2 


—7.56 


—3.33 


0.28 










CV Cha 


—11.17 


—3.18 


0.41 


—6.81 








CW Tau 


—4.81 


—4.17 


—0.49 


—7.80 


-1.94 


-7.24 


2 


DF Tau 


—3.97 


—4.33 


—0.63 


—6.93 


—1.08 


—6.49 


2 


DG Tau 


—7.16 


-3.82 


-0.17 


-7.49 


—11.52 


—19.14 


2 


DL Tau 


—12.20 


—4.15 


—0.47 


—7.41 


—14.91 


—24.63 


2 


DM Tau 


—20.73 


-5.20 


-1.42 


-8.48 








DoAr 21 


—4.62 


—4.49 


—0.77 


—7.78 








DoAr 24 E S 


-1.51 


-4.65 


-0.92 


-7.44 








DoAr 24 E N 


—4.02 


—4.52 


—0.80 










DoAr 44 


—3.68 


—4.92 


— 1.16 


—8.43 








DO Tau 


—4.18 


-4.37 


-0.67 


-7.28 


—2.28 


—9.56 


2 


DP Tau 


—5.78 


-4.81 


-1.06 


-8.24 


—3.40 


—7.50 


5 


DR Tau 


—6.86 


—3.90 


—0.24 


—6.82 


-7.39 


—16.78 


2 


DS Tau 


—6.99 


—4.75 


—1.01 


—8.22 


—1.23 


—4.96 


2 


EC 90 N 


—7.71 


—3.23 


0.37 










EC 90 S 


—2.21 


—2.87 


0.69 










Elias 29 


—3.10 


—3.29 


0.31 










EX Lup 


—3.70 


—3.83 


—0.17 










FS CmA 


—5.59 


— 1.66 


1.79 










FT Tau 


—6.91 


—4.86 


— 1.11 










FZ Tau 


—13.19 


-3.94 


-0.27 


-7.33 








CM Aur 


—29.90 


-4.51 


-0.80 


-7.66 


—7.12 


-12.20 




GQ Lup 


—3.38 


-4.85 


-1.10 


-8.15 








GV Tau S 


—8.32 


—3.49 


0.13 










GY 224 


—4.54 


—4.58 


—0.86 


-7.74 




—17.90 




Haro 1-1 


—46.90 


—4.50 


—0.78 










Haro 1-4 


—15.40 


—4.56 


—0.84 










Haro 6-13 


—6.70 


—4.11 


—0.44 


—7.02 


—4.49 


—5.61 


2 


HD 135344 B 


—1.99 


—4.86 


—1.11 


—8.35 


— 1.30 




3 


HD 141569 A 


—27.00 


-4.34 


—0.64 


—8.13 


—4.50 




3 


HD 142527 


—2.36 


-3.56 


0.07 


-7.02 


— 1.20 




3 


HD 142666 


—5.40 


—3.94 


—0.28 


—7.60 


—2.40 




3 


HD 143006 


—2.27 


—4.57 


—0.85 










HD 144432 S 


—5.52 


—4.03 


—0.35 


—7.69 


—5.10 




3 


HD 144965 


—24.96 


-3.85 


-0.20 










HD 150193 


-7.15 


-3.22 


0.38 


-6.97 


—5.50 




3 


HD 163296 


—7.68 


—3.33 


0.28 


—7.13 


—6.90 




3 


HD 169142 


—4.72 


—4.36 


—0.66 


—8.04 


—9.70 




3 


HD 179218 


—7.02 


—3.18 


0.42 


—6.76 


—10.10 




3 


HD 190073 


—5.52 


—2.11 


1.38 


—5.82 








HD 244604 


—5.28 


—3.48 


0.14 


—7.19 


—2.90 




1 


HD 35187 


—2.63 


—4.42 


—0.71 


—8.15 


-2.40 




1 


HD 37258 


—4.72 


—3.50 


0.13 










HD 37357 


—6.83 


—3.98 


—0.32 


—7.88 








HD 37411 


—4.10 


—3.42 


0.19 










HD 37806 


—5.19 


-2.59 


0.95 


-6.31 


—3.90 




1 


HD 38120 


—9.18 


—3.27 


0.33 


—7.04 


—9.60 




1 


HD 50138 


—6.96 


-1.70 


1.76 




-9.60 




1 


HD 58647 


—2.34 


-2.68 


0.86 


-6.36 








HH 100 


—3.71 


—3.40 


0.21 


—6.23 








HL Tau 


—6.00 


—3.74 


—0.10 


—7.01 


—4.71 


— 12.89 




IM Lup 


—5.63 


—4.82 


— 1.08 










IRAS 03259+3105 


-12.48 


-3.43 


0.19 










IRAS 03301-1-3111 


—7.07 


—3.95 


—0.29 










IRAS 034454-3242 


-0.51 


-4.18 


-0.50 










IRAS 042394-2436 


-9.56 


-3.87 


-0.21 










IRS 43 


-2.28 


-4.38 


-0.68 










IRS 44 E 


-2.20 


-4.55 


-0.83 










IRS 48 


-1.20 


-4.62 


-0.89 


-8.40 








IRS 51 


-5.53 


-3.95 


-0.28 










IRS 60 S 


-4.04 


-4.65 


-0.92 










IRS 63 


-6.26 


-4.24 


-0.55 










LkHa 25 


-9.49 


-3.11 


0.48 










LkHa 198 


-3.72 


-2.58 


0.95 










LkHa 220 


-10.80 


-2.86 


0.70 
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Star 


EWpf« 


loe: LT3f rt^ 


log -Lace 


log M^^ 


EW (Br7) 


EW (Pa/3) 


Rcfs 




fAl 
[^J 






[iU0 yr J 








LkHa 271 


-14.10 


-4.46 


-0.75 










LkHa 326 


-9.48 


-4.13 


-0.45 










LkHa 327 NE 


-6.41 


-4.09 


-0.42 










LkHa 327 SW 


-5.88 


-4.34 


-0.64 










LkHa 330 


-3.08 


-4.14 


-0.46 


-7.66 








LkHa 348 


-8.58 


-2.80 


0.75 










LLN 8 


-2.54 


-2.85 


0.71 


-7.14 








LLN 17 


-7.01 


-2.15 


1.34 










MCW 442 


-4.84 


-2.03 


1.46 










MWC 480 


-8.84 


-3.51 


0.11 


-7.21 


-7.00 




1 


MWC 758 


-3.59 


-3.74 


-0.10 


-7.35 


-5.40 




1 


MWC 1080 


-4.58 


-0.63 


2.72 










RNO 19 


-9.87 


-3.97 


-0.30 










RNO 90 


-8.07 


-3.85 


-0.19 


-7.40 








RNO 91 


-4.52 


-4.39 


-0.68 










ROX 43A 


-5.21 


-4.53 


-0.81 










RR Tau 


-5.18 


-1.96 


1.52 


-5.68 


-3.50 






RU Lup 


-6.91 


-4.15 


-0.47 


-7.75 








RXJ1842.9-3532 


-9.28 


-4.61 


-0.89 










RY Tau 


-4.06 


-3.56 


0.07 


-7.17 


-0.51 


-2.24 




S CrA A 


-6.06 


-3.77 


-0.13 


-7.36 








S CrA B 


-4.75 


-4.35 


-0.65 


-7.53 








SR 4 


-9.78 


-4.52 


-0.80 


-7.55 








SU Aur 


-6.16 


-3.74 


-0.10 


-7.31 


-0.67 


-0.63 




Sz 68 


-2.06 


-4.42 


-0.71 


-7.78 








T CrA 


-1.86 


-4.68 


-0.94 










T Tau S 


-6.34 


-3.03 


0.55 










T Tau N 


-5.51 


-3.58 


0.05 


-7.16 








TW Hya 


-27.38 


-5.07 


-1.30 


-8.67 


-7.27 






TY CrA 


-21.16 


-3.90 


-0.24 










UY Aur A 


-1.47 


-4.88 


-1.13 










UY Aur B 


-4.69 


-4.64 


-0.92 










UX Ori 


-4.32 


-3.32 


0.29 


-7.03 


-3.90 




3 


VSSG 1 


-3.35 


-4.60 


-0.87 


-7.56 




-8.90 


4 


VV CrA S 


-5.71 


-3.47 


0.15 










VV CrA N 


-2.46 


-4.36 


-0.66 










VV Ser 


-5.89 


-3.27 


0.33 


-7.13 


-13.00 






VW Cha 


-9.18 


-4.48 


-0.77 


-7.50 








VZ Cha 


-4.58 


-4.44 


-0.73 


-7.91 








V1685 Cyg 


-6.80 


-1.74 


1.72 










V1977 Cyg 


-3.00 


-3.11 


0.48 


-6.86 








Wa Oph 6 


-5.36 


-4.44 


-0.73 










WL 12 


-1.18 


-5.05 


-1.28 










WL 16 


-17.31 


-3.54 


0.08 


-6.77 




-8.40 


4 


WL 22 


-46.67 


-2.05 


1.43 










WW Vul 


-3.32 


-3.12 


0.47 


-6.86 


-5.40 




3 



References. — (1) Donehew & Brittain 2011 (2) Folha et al. 2001 (3) Garcia Lopez et al. 2006 (4) MuzcroUe et al. 1998 (5) Natta et 
al. 2006 

^ Adjusted for photometry and extinction, as described in text. Assumes distances in Table|3] If both NIRSPEC and CRIRES data were 
available, the average value is shown. 

^ Calculated assuming Lace = 0.8 GA/^- A///?* I IGuUbring et al]|1998l ). If not shown in Table [S] iJ* is derived from and T-,^ assuming 
L* = 4:-KR^aT^. Missing data indicate that either the stellar radius or mass is not available. 
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Table 5 

Literature Accretion Luminosities and Mass Accretion Rates 
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References. — (1) Andrews & Williams 2007 (2) Andrews et al. 2011 (3) Aspin et al. 2010 (4) Blondel et al. 2006 (5) Calvet et al. 
2004 (6) Donehew & Brittain 2011 (7) Eisner et al. 2005 (8) Fernandez et al. 1995 (9) Garcia Lopez et al. 2006 (10) GuUbring et al. 1998 
(11) Hartmann et al. 1998 (12) Hcrczeg et al. 2008 (13) Muzerolle et al. 1998 (14) Mcndiguti'a et al. 2011 (15) Natta et al. 2006 (16) 
Nisini et al. 2005 (17) Prato et al. 2009 (18) Valenti et al. 1993 (19) White & Chez 2001 (20) White & Hillenbrand 2004 

^ Description of abbreviations are as follows. UV Exc. = measurement of UV excess emission spectrum; Br 7 = use of Br 7-accretion 
luminosity relationship; Pa /3 = use of Pa /3-accretion luminosity relationship; U-Band Phot. = use of optical U-band photometric 
data; Balmer depth = based on relationship between the mass accretion rate and the excess in the Balmer discontinuity derived by 
IMuzeroUe eFall l l200i) 

^ This and all other Taurus accretion luminosities from I Valenti et al.l 11993) scaled to 140 pc distance. 
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